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TECHNICAL MEMORANDUM X-53148 


A FORTRAN PROGRAM TO CALCULATE AN ENGINEERING ESTIMATE OF THE 
THERMAL RADIATION TO THE BASE OF A MULTI-ENGINE 
SPACE VEHICLE AT HIGH ALTITUDES 


SUMMARY 


A method of estimating the radiant heat flux in a base of arbitrary 
shape from intersection regions caused by the interaction of hydrogen- 
oxygen engine exhaust jets is presented. An approximate method of 
generating the intersection region shape and tempera ture-pressufe pro- 
files is discussed. A computer program incorporating both of the above 
is described and instructions are given for its loading and use. The 
a accuracy of this program is expected to yield within an order of magni- 
tude of the true value of thermal radiation. 

I. INTRODUCTION 

A contributing factor to the base heating rate on upper stage 
vehicles is radiation from the exhaust plumes. Limiting the discussion 
to upper stage vehicles makes it unnecessary to consider the afterburning 
problem. With this simplification, this radiant heating problem lends 
itself to division into three distinct radiating regions: the super- 

sonic cores, the first order plume impingement regions, and the extremely 
complex flow field downstream from the base where the supersonic cores, 
first order regions, and higher order regions merge. The term ff first 
order plume impingement region” is used to describe the region created 
by the shock formed when two exhaust plumes intersect. This region is 
an area of high temperature and pressure. A second order region is 
formed by the intersection of two first order regions. This terminology 
can be extended to describe any number of intersecting regions and will 
be used throughout this report. 

Of the three types of radiating regions, the first two yield to 
analysis while the third, being virtually undefinable, does not. It was 
decided early in the analysis to consider the supersonic cores and first 
order regions as independent primary sources of radiation. To take 
account of the complex downstream flow field, it was decided to estimate 
a “lump sum” addition to the final overall radiative flux, this estimate 
being based on the numbers obtained for the primary radiation sources. 
Further analysis has revealed that of the two remaining regions, the 
supersonic cores are often insignificant; thus we have limited ourselves 
in the computer program described here to a calculation of the radiation 
from the first order intersection regions only. 



The supersonic cores in general are ac a considerably lover calcu- 
lated temperature than the intersection region and, in addition, are 
located closer to the base. Because of this, radiation is less intense 
and more completely blocked by components in the base region. However, 
fepr the user ' who feels that the supersonic cores will contribute 
materially to his particular problem, a program has been developed to 
calculate the additional radiation from these regions using methods 
identical to those described in this report. 

Proper use of the program discussed here will thus provide the user 
with an estimate of the incident radiation in a minimum of time and will 
show if a more detailed analysis is needed. 


II. CAPABILITIES AND LIMITATIONS 
A. Capabilities 


The program will estimate the radiation incident on a base of 
arbitrary configuration from any number of first order regions and will 
predict how much of this radiation will be blocked by intervening struc- 
tures in the base region. The user has the option of choosing the 
method of calculation of the temperature and pressure fields in the 
intersection region from the three methods programmed. In addition, 
the actual radiation is calculated by two different methods to allow 
comparison of two theoretical approaches to the problem. Of these 
methods, the average temperature and pressure method is more conservative 
since it assumes that the radiation passes through a transparent space 
from the radiating region to the base, when in fact it is passing through 
an absorbing medium. This method thus provides a probable upper limit 
to the heat flux rate. The other approach accounts for the absorption 
by an approximate method outlined by Hottel in Reference 1. Each method 
is presented with and without* a blockage correction factor. 

B. Limitations and Assumptions 

As mentioned earlier, this program computes only the radiation 
contribution from the first order intersection regions. These regions 
are assumed to have approximately elliptical cross sections and their 
temperature and pressure profiles are radially invariant. That is, 
temperature and pressure are assumed to vary in an axial direction only. 
These profiles and the geometry of the region are generated from input 
taken from a method of characteristics printout for a free jet expansion 
into quiescent air. Tne flow model for the first order regions is 
generated by assuming that this expanding plume impinges on a flat plate 
placed in the flow at a distance equal to one half the distance between 
nozzle centerlines. Oblique shock theory is then used to produce the 
geometry as described in Section IIIA of this report. 
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Both radiation methods use the emissivity data of Hottel given 
in Reference 1 without the correction for partial pressure. . These data 
are used in the form of a polynomial curve fit with a maximum error in 
fit of 3 percent. Radiation from each intersection region is assumed to 
pass through a transparent medium until it impinges on the base or is 
blocked by one of the solid components in the base region (such as 
motors and heat shield). 

Therefore, the program is restricted to use in cases where the 
only significant radiating species in the exhaust is water vapor. If 
the curves were fit to new emissivity data, the program could be changed 
to treat other radiating species. 

III. THEORETICAL BACKGROUND 


The gas radiation is treated with a comparatively crude approxi- 
mation. The emission and absorption coefficients, which actually vary 
in more or less unknown fashion through the emission spectrum of the gas 
are replaced by a lumped emissivity representing an integrated effect 
over the whole emission spectrum. This lumped or total emissivity was 
obtained from experiments on comparatively small samples of a homogeneous 
gas and extrapolated versus sample size, pressure, and temperature. The 
extrapolation becomes progressively more uncertain the further one moves 
away from the actual test conditions. Use of these data for estimating 
radiative heating of full scale rocket configurations should normally 
involve extrapolation of these data to larger beam lengths and lower 
pressures than those used in the basic experiments. Further, these 
measurements from homogeneous gas samples are used to estimate radiative 
emission from nonhomogeneous gas volumes. This may produce additional 
uncertainties . 

Two methods are proposed for this estimate: the average tempera- 

ture method and the absorption method. In both cases, the radiating gas 
volume is subdivided into conical pencil-like volume elements, extending 
along lines of sight from the base surface element in question. Tempera- 
ture, density, and partial pressure of the radiating species vary along 
the gas pencil. 

For the average temperature method, the properties of the gas, i.e., 
T 4 and partial pressure, are averaged over the length of the pencil. 

The "pencil” is then treated as a homogeneous gas sample, yielding a 
corresponding emissivity with a corresponding view factor. From these, 
the energy flux contribution that impinges on the base surface element 
under investigation is computed. 



In the absorption method, the "gas pencil" is divided longitudinally 
into individual elements of different pressures and temperatures. The 
energy radiated from each longitudinal element toward the base surface 
element under investigation is subject to absorption losses in the 
pencil elements it has to pass through. This absorption loss is esti- 
mated by the assumption that the whole gas pencil between the radiating 
element and the surface of the intersection region has the same gas 
properties (T, Ph 2 o) as the radiating element. While this model admittedly 
loses any effects of property variation on the absorption, it significantly 
reduces the complexity of the estimate by eliminating one integration. 

To simplify the presentation, the following discussion will consider 
the program as though it were divided into four separate programs. 


A. Geometry and Temperature-Pressure Profiles 


From a method of characteristics analysis of a free plume expan- 
sion into quiescent air, one may obtain values for Mach number (M) and 
flow direction angle ( 0 ) at various points in the exhaust plume. Values 
of M and 0 along a line of constant radius are interpolated from this 
output (Figure 1). In addition, one obtains the plume expansion radius 
(radial distance from nozzle centerline to jet boundary streamline) as 
a function of axial distance. Using these input data, the initial region 
geometry and temperature-pressure profiles are generated using the follow- 
ing procedure. 

1 . Geometry of the Intersection Region 

In Figure 2, a representative cross section is drawn. The 
semi-major axis of the cross section, ADIS, was determined by connecting 
the overlap of circles whose diameters were established by the free plume 
boundaries from the method of characteristics. Determination of the semi- 
minor axis was more difficult. It was assumed that the flow field from 
two adjacent streams deflects by the same amount. This assumption pre- 
serves symmetry about the major axis of the region and is compatible with 
a hypothetical model consisting of a plane placed half v;ay between two 
adjacent engines. This model is sketched in Figure 3. At points 0, 1, 

2, 3, etc., the Mach number and flow angles are known. The flow imping- 
ing on the plate at point 0 produces an oblique shock at an angle S 0 
given by the equation 
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where M x , the Mach number before the shock, is the Mach number M 0 which 
we have obtained from the method of characteristics analysis. This shock 
is shown as OB-1. It is clear that it is at an angle of (5 - 0) from the 
hypothetical flat plate. The mass flow between streamline SL-0 and SL-1 
is now assumed to be turned at the shock so that it is flowing downstream 
parallel to the wall. Because of continuity, the flow occupies a certain 
distance from the plate to some streamline SL-1 1 . Next, the Mach number 
and flow angle known to exist at point 1 are transferred radially to 
point l 1 . This transfer is justified on the grounds that the temperature 
and pressure are assumed constant in a radial direction and that thus the 
flow angle and Mach number must be constant also. At point l 1 a new 
shock angle S x is calculated as before awd another section of the inter- 
section region boundary is determined. The process is repeated until 
the entire region geometry is obtained. It is immediately obvious that 
this approach ignores mass flowing between point A and point l 1 . This 
assumption should be reasonable if the value of (5 - 0) o is small and 
the value of (0) x is large. Tnis case corresponds to the highly expand- 
ing plume, which is the case normally encountered in high altitude 
operation. 


Equation (1) is limited in that it can only generate slock 
angles for values of flow angle less than a certain value. Figure 12 
shows a parametric plot of this equation for a j value of 1.23 and for 
Mach numbers ranging from 1.1 to 20. It is clear that since plume Mach 
numbers seldom exceed 10, flow angles of 55 degrees and above will not 
yield solutions to equation (1). In these cases, a shock angle of 
65 degrees was assumed since this is approximately the maximum shock 
angle that can be produced by equation (1) for Mach numbers in the range 
of those encountered in jet plumes. (Note that there are two solutions 
to equation (1) for any flow angle. The highest, or "strong shock" 
solution, is not used in the analysis.) 

The process outlined above is repeated until a complete 
geometrical profile of the region is developed. Figure 4 illustrates 
the application of the procedure for several increments. Referring to 
the Figure, we can see that the geometrical relation which links any one 
radial distance to the preceding radial distance is 


(B/D) n = (B/D) n _ i + tan (5 - 0) n _ 1 (X/D) n - (X/D) n _ r (2) 
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2. Pressure-Temperature Profiles 


With the geometry of the intersection region computed, the 
program returns to the input data. The pressure- tempera ture profiles 
are generated using the same model as before. Postulating as before 
the information of an oblique shock whose angle is determined by the 
input flow angle for that particular axial distance, the proglem resolves 
itself into the determination of static temperatures and pressures behind 
a shock from the total conditions before the shock. Three methods of 
calculating these temperatures and pressures are programmed and may be 
selected by the user. In general, the oblique shock temperature and 
pressure equations should be used since they are more conservative than 
the other two models and because their use maintains consistency with 
the equations used to generate the region geometry. 

Using the subscripts 1, 2, and 0 to denote, respectively, 
static conditions before the shock, static conditions after the shock, 
and total conditions, the following isentropic perfect gas relation 
holds: 



Also the adiabatic perfect gas relation is 



The oblique shock relations used are taken from NAGA Report 1135 as 


Is = [2y Mi 2 sin 2 5 - (r-- l)][(y - 1) M n 2 sin 2 8 + 2] . 

Ti (y + Vp HjZ sin k; 6 W 


and 


P p _ 2y M-) 2 sin 2 & - (y - 1) 
Pi " 7+1 


6 



As mentioned in the geometry section, cases exist where the 
shock angle 5 cannot be determined. In these situations, the program 
automatically switches to another method of calculation using the Rankine- 
Huginot equation 


To - T, 


L(1 






Mi 2 ) (l + a) + p + Mi 2 


(7a) 


where 


a 


= 7 1 

y + 1 


and 


P = 7 M-l 2 sin 2 e 


(7b) 


for the temperature calculation. 

For the pressure calculation, the modified Newtonian Impact 
Theory equation 


f 2 = C (//2) M-j 2 + 1 
*1 P 


( 8 ) 


and 


C = 2 sin 2 0 (9) 

P 


is used. These equations may also be called for by the user if desired 
by inputing a control value as described in Section VIIIA of this report. 


In addition. 


the isentropic flow equation 



X 

2 


M x 2 + 1 



( 10 ) 
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is programmed and may be requested by use of the control variable. The 
pressure calculation is made using equation (8) . 

These calculations are repeated progressing -in an axial 
direction until the entire profile has been generated. Due to the 

assumption that the static temperature before the shock is constant, 

it follows that the static temperature behind the shock does not vary 

from point 0 to point 1* in Figure 3 and thus also from 0 to 1. Then, 

at point 1, a step discontinuity will exist and the temperature from 1 
to 2 will be the temperature calculated from the flow angle at point 1. 
If a plot of temperature vs axial distance were obtained, these discon- 
tinuities would cause the curve to resemble a "step" function. To 
smooth the curve, linear temperature gradients and not constant tempera- 
tures are assumed to exist between analyzed axial points. 

B. The Average Temperature Method 

The Average Temperature Method is an analysis of the problem 
which uses the fundamental Stephan-Boltzman Law 


Q = aT 4 . 


( 11 ) 


Obviously, a sort of mean temperature and pressure for each radiating 
increment must be developed. The method of calculating these mean values 
is described below. 

The intersection region was first divided into a number of 
increments. A line of sight from a point on the base to the lower edge 
of each increment was established and extended through to the rear face 
of the plume. A view factor, F, to account for the fact that the point 
on the base only "sees’ 1 one side of each radiating increment, was calcu- 
lated following the analysis of Sparrow [2]. This derivation is per- 
formed in detail in Appendix A. This factor is calculated for each 
individual increment radiating to a point (Figure 5). 

The mean temperature along the beam length itself is 



(T (x) ) 4 dx 



( 12 ) 
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where x x and x 2 are the axial coordinates of the points where the beam 
enters and leaves the intersection region, respectively. A similar 
formula is applied for the evaluation of the average pressure along 
the beam length. Using this average pressure, the partial pressure of 
the water vapor in the exhaust is 


V ■ (FX1) P AV> 


( 13 ) 


where FXl is the mole fraction of H 2 0 vapor in the exhaust plume. 

With the partial pressure, temperature, and beam length deter- 
mined, the emissivity is obtained from the data of Hottel [1]. These 
emissivity data have been curve-fitted with a series of polynomial 
expressions such that the maximum error is ±3 percent in fit. From this 
the emissivity may be determined and used in the final expression for Q, 
the heat flux 


Q = oF e T 4 (14) 


C. The Absorption Method 

The absorption method employs the equation developed by Hottel 
[1] to account for the absorption in the intersection region. The equa- 
tion employed is 


Q 



go X 


a T 4 (de/dX) cos 0dto dX. 


(15) 


In this expression, X is the optical beam length and is evaluated as 



(16) 
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and 


cos 0 dco ( 17 ) 

Cd 

is a solid angle shape factor previously called F. Then 

q ' fo / t4 a<ft) d(pL) - < is ) 

P„L 

The plume is divided as before into increments and a pencil 
ray (line of sight) erected to the lower edge of each increment. This 
line of sight is now the optical beam along which equation (18) must 
be evaluated. To do this, we must first evaluate the pressure and 
temperature at various points along the beam. This is done using a 
straight line interpolation method which is quite valid because of the 
smooth nature of the temperature and pressure profiles and to the small 
distance between points at which the evaluation is made. 

The emissivity slope is easily evaluated since the curve fits 
mentioned in Section IIIB are functions of (PL) with T as a parameter. 
Thus, an explicit differentiation of these polynomial curve fits is 
possible without resorting to numerical techniques. Since all the 
variables are functions of the position variable L, we may now proceed 
to simplify equation (18). 

PL 2 • 

Q = Fa J T x 4 

PLi 

s ince 


d (PL) = PdL (20) 


(l/ p ) |f PdL, 


(19) 
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by the assumption that temperature and pressure are constant in the 
particular increment of length under analysis. Then 


Q = Fa 



PLi 


( 21 ) 


Equation (21) is evaluated using Simpson’s rule and the value for F 
already computed in Section IIIB. 

D. Blockage 

In certain base configurations, a great deal of the radiation 
from the exhaust plumes will impinge upon the motors, heat shield, and 
other obstructions in the base, and will not therefore contribute to 
the heat flux rate to the base. To evaluate how much of the radiation 
is blocked, a section of the program generates a blockage factor and 
applies it to each intersection region. This blockage factor is not 
related to the form factor F computed in Section IIIB. The form factor 
expresses what portion of the intersection region the point "sees" 
assuming a transparent space between radiating region and point. The 
blockage factor increases the realism of the model by removing the 
"transparent space" assumption and substituting the actual base obstruc- 
tions. 


It was decided that the easiest method of generating a blockage 
factor was to determine if each individual beam was blocked and if it 
were to consider the radiation from the increment associated with that 
beam to be 0. Since the increments have widths which in some cases are 
appreciable, three lines of sight are erected to each increment, one to 
each side and one to the center. Each line is assumed to represent one- 
third of the radiation incident from the region. Thus, blockage of each 
beam reduces the incident radiation by one-third. 

The obstructions are first represented as circles in a three- 
dimensional, Cartesian coordinate system located as follows. 

1. The xy plane is located at the exit plane of the motors 
such that the z coordinates of all obstructions and the points to be 
analyzed are negative or zero. 
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2. The xz plane passes through the center of the heat shield 
and through the point being analyzed with the negative axis in the direc- 
tion of the point being analyzed such that "x" and "z" coordinates of 
the point are negative and the ,f y ,f coordinate is always 0. 

3. The yz plane is oriented as usual in a right-handed Cartesian 
coordinate system. See Figure 6 for a graphical illustration of this 

sys tem. 


The obstructions are assumed to be circles located in, or 
parallel to, the xy plane. In this way, a three-dimensional nozzle can 
be fairly well represented by, for example, three circles with radii 
equaling the radii of the nozzle at three different heights. The model 
expressing this nozzle would then become three circular areas in space 
given by the relations (in our coordinate system) 


(x - h) 2 + (y - k) 2 g r z 2 

(x - h) 2 + (y - k) 2 ^ r 2 2 

(x - h) 2 + (y - k) 2 g r 3 2 


z 


z i 


z - z 2 


z = z 3 . 


( 22 ) 

(23) 

(24) 


If the obstruction has its axis parallel to the axis of the vehicle 
(which is our z-axis) then ,f h H and ,f k M are constant for the three equa- 
tions. The "h" and ,, k tl coordinates are the M x M and M y M coordinates of 
the center of the obstruction. If the axis of the body is tilted with 
respect to the vehicle axis then "h" and ,! k M vary. They are then the 
coordinates of the center of the obstruction at that particular height. 

A slight error is thus introduced because of the assumption that the 
circles lie parallel to the xy plane when, in fact, they are tilted to 
be perpendicular to the obstruction axis; but this error will not be 
appreciable unless the tilt is considerable. 

With equations (22), (23), and (24) defining three areas in 
space, it is clear that, for a line of sight to avoid intersecting one 
of these areas, yet in fact to be blocked by the nozzle, the line of 
sight must be almost parallel to the xy plane. Since in practice the 
intersection region must begin at some finite distance above the nozzle 
exit plane, and the point to be analyzed must be some distance below the 
nozzle exit plane, the line of sight angle should be sufficient to allow 
these three circles to adequately describe the blockage. In many cases, 
in fact, one circle suffices. (In theory, a large number of circles 
used as above could completely describe an axisymmetric obstruction. 

In practice, however, only forty equations may be used in the program 
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to describe the whole base or a portion of it if so desired. Therefore, 
one should not use more circles than are necessary on any one obstruc- 
tion.) 


Once the entire base picture has been presented in the 
simplified terms shown above, the coordinate system is shifted so that 
it is now centered at the point to be analyzed, the orientation of the 
axes being as before. Now it is clear that all the obstructions are at 
positive "z" coordinates and most are also at positive "x" coordinates. 

To illustrate this transformation, assume a point P located at 
(x' , 0, z' ) as shown in Figure 6 where z = 0 is defined to exist at the 
nozzle exit plane. Assume also a point in an intersection region located 
at (x^, y^). Its z coordinate is unimportant. Assume also a circular 
blockage located at (h, k, z Q ) with a radius "r." The equations for 
this blockage region are 


(x - h 2 ) + (y - k) 2 g r 2 . z = z 0 


(25) 


To transfer these equations to a coordinate system centered at P, the 
following transformation equations are needed. (The subscript 2 refers 
to the new system.) 

x 2 = x - x' y 2 = y z 2 = z - z' . (26) 


Then, the transformed blockage equation is 


(x 2 + x' - h) 2 + (y - k) 2 g r 2 z 2 = z Q - z' . (27) 


Given a line of sight with an angle, p, from the horizontal 
(the new xy plane) and knowing the z coordinate of the obstruction (Fig- 
ure 7), we may calculate the distance (HYP) from P to the point where the 
line of sight passes through the plane of the obstruction (z 2 = z Q - z'). 


HYP 


z ? 

tan (p) 


(28) 
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Also, from the coordinates of the center of the intersection region 
which must be crossed by the beam, an angle 0 may be calculated as 
follows: 


0 = tan 



( 29 ) 


With 0 and HYP, the coordinates of the point where the line of sight 
passes through the intersection region may be calculated. 

x Q = (cos 0) HYP (30) 

y o = (sin 0) HYP. (31) 

The sign of y Q is the same as the sign of y^ since P is on the x-axis. 

Now, if tue calculated values of x Q and y Q satisfy equation (25) 
the line of sight is blocked. If not, the process is repeated with new 
parameters. Finally, a new summation of all the unblocked radiation is 
made to give a better estimate of the environment. 


IV. PROGRAM DESCRIPTION 

The program is written in Chain Fortran for the GE 225 computer. 

The use of Chains was caused by the large number of subscripted variables 
and makes the program a compendium of four smaller programs. Each of 
these smaller programs contains the coding of the equations contained in 
one of the four sections just discussed. 

This segmenting makes checkout and error investigation quite simple 
since one can immediately isolate a given section as the one in which the 
error occurred and restrict his investigation to it. 

Each of the four sections follows closely the pattern laid down in 
the preceding portion of this report. Numerous interpolation, approxi- 
mation, and averaging methods which are similar to techniques developed 
in the standard texts [3] are used as needed in the program. No 
unusual numerical differentiation and integration techniques are 
employed, and nothing other than simple geometry is used in the various 
form factor and beam length calculations. In short, the theory pre- 
sented in the preceding portion of this report is followed closely with 
little modification. 
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Figure 8 contains a printout of the program as it is read into 
the machine. Notice the marginal notations separating the four sec- 
tions of the program. If each of these sections is compared with the 
corresponding theory section of this, report by a person familiar with 
Fortran II, there is no reason why he should not be able to follow the 

programming and to verify his results in any way he wishes. 

Figure 9 contains a flow diagram to completely illustrate the flow 

, of control throughout the program. 


V, NOTATION 


Program Symbol Description and Units 

(J) Numerical subscript where J has a value 

(or a similar subscript) from one number to another (i.e., 1 to L) 


A (I, J) 

(J = 1 to 4) 


Table of input from the method of 
characteristics . 


A (I, 1) 


A(I, 2) 


CHAIN I 


A (I, 3) 


A(I, 4) 


Distance, nondimens ionalized by DIA, 
measured along the nozzle axial line 
from the exit plane of that nozzle. 

The 'first value must be the initial 
intersection point. 

Mach number corresponding to A(I, 1) 
but located on the axial line of the 
intersecting region (not to be confused 
with critical condition Mach number, 

M*). 

Flow direction angle corresponding to 
A (I, 1) but located on the axial line 
of the intersecting region, in radians. 

Radius, nondimens ionalized by DIA, of 
the circle formed by a plane drawn perpen 
dicular to the nozzle axial line and 
cutting the exhaust plume for a 
corresponding A(I, 1). 
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Description and Units 


Program i 

\ 

A (I, 1)’ 
A (I, 2) 
A(I, 3) > 
A (I, 4) 

ABSIS 

ADIS ( J) 

ANGLE (J) 
AVTEM(J) 
AVP(J) 
BDIS (J) 

BLED (J) 

C 

DELI 

DELM 


nnbol 


CHAIN IV 


Center and radii of the blockage regions 
located in the (x, y, z) coordinate system 
illustrated in Figure 6, ft. 


An increment length, nondimens ionalized by 
DIA, that is 1/50 that of the beam length, 
FLE(J). 

One-half the length, nondimens ionalized by 
DIA, of the major axis of the elliptical 
plane for the intersecting region. 

The angle formed by the beam and its pro- 
jection in the base plane, radians. 

Average temperature for the beam length, 

FLE(J) , °R. 

Average pressure for the beam length, FLE(J) , 
AIM. 

One-half the length, nondimens ionalized 
by DIA, of the minor axis of the elliptical 
plane for the intersecting region. 

Dimensional value for the beam length, ft. 

The angular difference between the shock 
angle, DELTA, and a corresponding flow 
direction angle), radians. 

Minimum shock angle for a given GAMMA and 
Mach number, FM(J) (zero flow direction 
angle) , radians . 

Maximum shock angle for a given GAMMA and Mach 
number, FM(J), radians. 
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Program Symbol 


Description and Units 


DELTA Shock angle for a given GAMMA, Mach number, FM(J), 

and flow direction, T(J), radians. 

DIA Diameter of the nozzle exit, ft. 

EMIS Emissivity of the beam length, FLE(J). 


F 

FINTX 

FIFTY 

'FK 

FLAG 

FLE(J) 


Shape factor corresponding to a given beam. (F 
is the view factor associated with each beam 
length that radiates energy back to the viewpoint.) 

x-coordinate of intersection region center in 
coordinate system of Figure 6. 

y-coordinate of intersection region center in 
coordinate system of Figure 6. 

The distance, nondimens ionalized by DIA, in the 
base plane from the point of view to the axial 
line for the intersecting region. (The value of 
FK is always negative.) 

Value, either +1 or 0, instructing computer on 
how to proceed. 

Length, nondimens ionalized by DIA, of the beam 
within the bounds of the intersecting region. 


FM(J) Mach number on the axial line of the intersecting 

region corresponding to XD(J) (not to be confused 
with critical condition Mach number, M*) . 

FXI Mole fraction of each of the gas components in the 

exhaust plume. 

GAMMA Ratio of specific heats for the exhaust components. 

L Number of desired data values for XD(J). 


LL 

LLLL 


Number of data values input for each A (I, 1), 
A (I, 2), A (I, 3), and A(I, 4) in Chain I. 

Number of blockage region circles. 
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Program Symbol 


Description and Units 


MCON 


PBPX 

PHI 


PH20 

PO 

P(J) and FW 

PX 

PY 

PZ 

QRAD 

QRADS 

QRDA 

QRDAS 

R 


If +1: oblique shock theory temperature and 

pressure used. 

0: Rankine-Huginot temperature and Newtonian 

impact theory pressure used* 

-1: isentropic temperature and Newtonian impact 

theory used* 

Ratio of base pressure to exit pressure. 

View angle formed by one-half the minor axis, 
BDIS(J), of the elliptical plane for the inter- 
secting region and a corresponding projection of 
part of the beam into the same elliptical plane 
(Figure 6), radians. 

Partial pressure for the beam length, FLE(J), ATM. 

Chamber pressure of the nozzle, ATM. 

Static pressure for a location in the exhaust 
plume, ATM. 

x-coordinate of point being analyzed in coordinate 
system of Figure 5. 

y-coordinate of point being analyzed. 

z-coordinate of point being analyzed. 

Radiant heat flux by the average temperature and 
pressure method for the beam length, FLE(J), 
BTU/ft 2 -sec. 

The total accumulated heat flux contribution by 
average temperature and pressure method for two 
or more beam lengths, BTU/ft 2 -sec. 

Radiant heat flux by absorption method for a 
single beam length, FLE(J), BTU/ft 2 -sec. 

Total accumulated radiant heat flux by absorption 
method for two or more beam lengths, BTU/ft 2 -sec. 

One-half the distance, nondimens ionalized by DIA, 
between two adjacent nozzle center lines. 
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Program Symbol 


Description and Units 


RC(J) 

RD(J) 

SABS IS 

SSF 

SXD 

TAB 

TABL 

TDIS ( J) 

TEMP(J) 

THETA 

T(J) 

TO 


One-half the length, nondimens ionalized by DIA, 
of an axis in the elliptical plane for the 
intersecting region that is formed by a beam where 
it pierces this region and a point on the axial 
line of the same region. 

Radius, nondimens ional ized by DIA, measured in 
the same manner as that of A(l, 4) and correspond- 
ing to XD(J). 

Total accumulated ABSIS increments to point in 
question. 

The total accumulated shape factor for two or 
more beam lengths. 

Perpendicular distance, nondimens ional ized by 
DIA, between the base plane and the nozzle exit 
plane. 

Table of flow direction angles, radians, obtained 
from a plot of flow direction angle, shock angle, 
and Mach number (Figure 13). 

Table of Mach number values obtained from a plot 
of flow direction angle, shock angle, and Mach 
number (Figure 13) . 

One-half the length, nondimens ional ized by DIA, 
of an axis perpendicular to RC(J) and in the same 
elliptical plane for the intersecting region. 

and TW Static temperature for a location in the exhaust 

plume, °R. 

Flow direction for a given GAMMA, Mach number, 
FM(J) , and shock angle, DELTA, radians. 

Flow direction angle on the axial line of the 
intersecting region corresponding to XD(J), 
radians . 

Chamber temperature of the nozzle, °R. 
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Program Symbol 
XD (J) 


XD (L + 1) 


XO 


Description and Units 

Arbitrary distance, nondimens ionalized by DIA, 
that is measured in the same manner as is A(I, 1). 
(Simply numbered values for XD(J) make it easier 
for the user to plot the results. XD(J) can be 
set equal to A(I, 1) if so desired.) 

Selected value, nondimens ionalized by DIA, that 
is slightly larger than the last input value for 
XD(J). (The value, XD(L +1), is necessary so 
that the solution of shape factor can be completed 
for the last increment.) 

x-coordinate of point at which line of sight 
to I.R. pierces plane of obstruction. 


YO 


y-coordinate of point at which line of sight to 
I.R. pierces plane of obstruction. 


ZO 


z-coordinate of point at which line of sight 
pierces plane of obstruction. 


VI. INPUT 


The program requires as input certain characteristics of the rocket 
engines on the missile in question, the base geometry and a flow field. 
Since the input required is quite extensive, it will be treated as three 
separate topics. In addition, the discussion will be restricted to the 
input and how it is obtained leaving the actual load format for a separate 
section. 

A. Flow Field Description 

The program requires extensive input data to describe the 
downstream flow field. This input is received in the form of tables 
of the free exhaust plume properties, Mach number * M, flow angle, 0, 
and plume radius, R, as functions of axial distance. The axial distance 
and the plume radius are nondimens ionalized by the exit diameter and 
will hereafter be referred to by their symbolic names X/D and R/D. The 
Mach number is not the characteristic Mach number M* but is related to 
it by 


2M* 2 

M n/ (7 + 1 (7 - 1) M*^ 


(32) 
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The angle q is expressed in radians. These tables are obtained from 
a method of characteristics analysis as follows. 

In Figure 1 it may be seen that by a two-way interpolation 
into the method of characteristics output, values of M and 0 may be 
found along a line of given radial distance from the nozzle centerline. 

In addition, at the specific X/D location at which the M and 0 1 s are 
found, the free plume radius must also be input as shown in the figure. 

The first X/D chosen must be at the initial intersection point and 
additional X/D*s down the plume may be chosen arbitrarily; however, 
the bottom number must not exceed 60. 

Given now this set of tables of M, 0, and R/D as functions of 
X/D, the user next inputs the X/D distances at which he wants his lines 
of sight to pierce the intersection region. The input fixes the number 
of increments into which the intersection region will be divided; thus, 
it will generally be advantageous to input as many as practical. The 
number of these points is in no way related to the number of points in 
the above tables except that no more than 60 may be input, and the 
highest and lowest value must be within the range of the previously 
input tables. 

A review of the oblique shock equations will show that the 
equation relating the shock angle 8 to the flow direction angle 0 does 
not allow a 8 to exist for all 0 values. In particular, for a given 
specific heat ratio and Mach number, there exists a maximum 0 beyond 
which a corresponding shock angle does not exist. The program is made 
aware of these limits by inputting a table of Mach numbers and the 
corresponding maximum 0 values so that it will not try to solve for a 
nonexistent angle. This table consists of ten values and must be 
obtained for the correct specific heat ratio as it varies strongly 
with that parameter. 

B. Rocket Motor Input 

The given rocket motor must be characterized for the program 
by inputting its total temperature, total pressure, its exit diameter, 
the specific heat ratio of its exhaust gases and the pressure ratio 
(exit pressure to ambient pressure) at which the method of characteristics 
analysis was run. 

C. Base Geometry 

Referring to Figure 10, the program needs the distances SXD 
and FK, both in the nondimens ionalized (divided by the exit diameter) 
mode. In addition, the location of the blockage regions and intersection 
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regions in a coordinate system oriented as described in Section HID are 
needed. The dimensions of this coordinate system must be in feet. 

In addition to the above input, a number of constants which 
effect control of the execution cycle are needed and are described later 
in the load format section. 


VII. OUTPUT 

The output is limited since execution time is already quite long 
and excessive print statements only make it longer. As may be seen in 
Figure 11, the output is arranged in the same logical order as was 
followed in Section III. The first section is devoted to the tempera- 
ture and pressure profiles; the second to one radiation method; the 
third to the absorption method; and the fourth to the blockage correction 
factors. As may be seen, the total cumulative radiation, as well as the 
radiation from the individual increments, is summed in each case. 

Note that in the blockage correction section each increment is 
divided into left, middle and right sides. The coordinates at which the 
line of sight to each of these three points pierces the planes of the 
various obstructions (x Q , y Q ) are printed out. If the line of sight is 
blocked, the last set of x Q , y Q coordinates printed out are within a 
blockage circle. This is shown on the output in Figure 11. 

Figure 11 contains the complete output for the sample problem in 
Section VIIIB. 


VIII. LOAD FORMAT AND SAMPLE PROBLEM 
A. Load Format 


A few basic rules for writing a load format for the GE 225 
Fortran are as follows: 

1. Fixed point numbers must have an X before the number and 
must be immediately followed by a comma (i.e., X31,). M X31 , ,f 
will be misinterpreted. 

2. Each value must be separated by a comma (i.e., .21, .2646, 

• • • ) . 

3. Every data card must end with an asterisk (i.e., 

.28, .29, .30*). 
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4. Floating point numbers can be written either in decimal form 
without exponent or with exponent vi.e., .000539 or 5.39E-4). 

5. Any one data cara cannot^ contain values listed on more than 
one read statement (RCD) . 

6. Data must be set up in the same order as it is listed in 
the read statement (RCD). 

7. All floating point numbers must have decimal points (i.e., 
the expression 31 is not permitted) . 

These rules must be followed explicitly. 

In the following discussion the input will be separated into blocks. 
Each block of data must go on one or more cards. The first item in each 
new block must b^ the first item on a new card; i.e., no card may contain 
data from two blocks. 

a. Block 1 

(1) A fixed point number expressing the number of values 
in the input table of M vs X/D. This number is of course also the 
number of values in the 0 vs X/D table and the R/D vs X/D table and must 
not exceed 60. 


(2) A fixed point number giving the number of X/D loca- 
tions at which the user desires to irtcrement his intersection region. 

(3) A fixed point number with its value determined from 
the following. 

(a) X + 1 if oblique shock relations are desired 
in the solution of the temperature and pres- 
sure profiles. 

(b) X + 0 if Rankine-Huginot relations are desired. 

(c) X - 1 if isentropic flow relations are needed. 

(4) A floating point number giving one-half the distance 
between the centerlines of the two motors which are generating the inter- 
section region. This number is the actual distance nondimens ionalized 
by the exit diameter of the motor. 

(5) The value of the chamber total pressure; atmospheres ~ 
floating point. 
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(6) The value of the chamber total temperature; °R ~ 
floating point. 

(7) The exhaust gas specific heat ratio; floating point, 

b. Block 2 


A floating point number expressing the base pressure to 
exit pressure ratio for which the method of characteristics run was made. 

c. Block 3 


A set of floating point numbers in ascending order express- 
ing tne X/D coordinates of the points at which the M and 0 will be given. 

d. Block 4 


A set of floating point numbers expressing the Mach numbers 
corresponding to the X/D locations given in Block 3. 

e. Block 5 


A set of floating point numbers giving the corresponding 0 
values; radians. 

f. Block 6 


A set of floating point numbers giving the corresponding 
R/D f s. (Note: Together Blocks 3, 4, 5, and 6 give the columns of a table 
into which the machine may go to interpolate the value of any parameter 
(M, 0, R/D) at any given X/D value during later execution stages. As 
a result, extreme care must be taken to assure that proper correspondence 
is maintained between the elements in the columns.) 

g. Block 7 


A series of floating point numbers giving the X/D locations 
at which the plume is desired to be incremented. The number of these 
points is given in Block 1, part 2. 

h. Block 8 and Block 9 


These two blocks contain the maximum values of M and 0, 
respectively, which are used to prevent the program from solving for a 
nonexistent shock angle. Referring to Figure 12, notice that for a 
Mach number 4.0 and a flow angle 0 of 55 degrees no shock angle 5 exists. 
Thus, to prevent the program from attempting to solve for shock angles 
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in cases like this, a table of Mach numbers and the maximum corresponding- 
0 angles is input* Figure 12 shows such a table far Figure 12 # This 
table must have ten values - 

i* Block 10 

(1) Exit diameter of the rocket engines used on the base 
being considered. A floating point number ~ feet. 

(2) The percent water vapor in the exhaust plumes. Float 
ing point ~ dimensionless. 

(3) A selected X/D location equal to the last X/D location 
in Block 7 plus .1. Floating point ~ nondimens ionalized. 

j. Block 11 

(1) The distance in the plane of the point being analyzed 
from the point of view to the axial line of the inter- 
secting region. This value is always given a negative 
sign. A floating point number ~ nondimens ionalized 
(see Figure 10) . 

(2) The angle formed by the semi-major axis, BDIS(J), of 
the elliptical plane for the intersection region and 
a corresponding projection of part of the beam into 
the same elliptical plane. A floating point number 
~ radians (see Figure 14) . 

(3) The perpendicular distance between the plane of the 
point being analyzed and the nozzle exit plane. A 
floating point number ~ nondimens ionalized (see Figure 10). 

k. Block 12 

(1) The number of blockage circles being considered. Must 
not exceed 40. A fixed point number ~ dimensionless . 

(2) The x-coordinate of the center of the intersection 
region being considered. A floating point number 
(including sign) ~ feet. 

(3) The y-coordinate of the center of the intersection 
region being considered. A floating point signed 
number feet. 
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1. Block 13 


A series of numbers expressing the x-coordinates of the 
blockage circle centers. Floating point signed numbers ~ feet. 

m. Block 14 

A series of numbers expressing the y- coordinate of the 
blockage circle centers. Floating point signed numbers ~ feet. 

n ♦ Block 15 

A series of numbers expressing the z-coordinate of the 
blockage circle centers.' Floating point signed number ~ feet. 

o. Block 16 

A series of numbers expressing the radius of each of the 
blockage circles. A floating point number ~ feet. 

p. Block 17 

(1) A value, either +1.0 or 0.0 depending upon whether a 
new intersection region is to be analyzed next or just 
a new point. This number is a command to the machine. 

+1.0 if a new point is to be analyzed or 
if this is the last piece of data. 

0.0 if a new intersection region is to be 
cons idered. 

(2) This value is +1.0 if the intersection region being 
analyzed has an x-coordinate greater than the x-coordin- 
ate of the point being considered. If the intersection 
region has an x-coordinate less than the point being 
considered, then the value in this position is the 
x-coordinate. A floating point negative (all x-coordin- 
ates are negative according to the defined coordinate 
system in Section HID) number ~ feet. As clarifica- 
tion, a case where the x-coordinate must be input in 
this manner is illustrated in Figure 15. 
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B. Sample Problem 


As a sample problem an analysis of the radiation incident on 
a point on the thrust structure of the S-II stage from two of the eight 
I.R.’s will be made. Figure 16 shows the base with some of the necessary 
geometrical input. Figure 17 shows the same base drawn in the necessary 
coordinate system to input the coordinates of the intersection regions 
and blockage regions. Other necessary input data concerning the engine 
and other features of the base are: 

Chamber Total Pressure: 43.004899 atmospheres. 

Chamber Total Temperature: 5760 °R 

Exit Diameter: 6.6666 feet 

Distance between Engine Centerlines: 12.37 feet 

Specific Heat Ratio of Exhaust Gases: 1.23 
Percent Water Vapor in Exhaust Gases: 63% = .63. 

Now a method of characteristics analysis of a J-2 engine free 
plume expansion into quiescent air must be obtained. From this we 
develop the following table. 


TABLE I 

(Note that the axial distance and radius have already been divided 
through by the exit diameter, 6.6666 ft.) 


A 

B 

c 

D 



1 

1 r 1 

X/D 

M 

Theta - 0 

Free Plume 

Expansion Radius - R/D 

55300 

5.55 

.59932 

.932 

60577 

5.678 

.61166 

.971 

68793 

6.121 

.69823 

1.022 

75513 

6.027 

.65766 

1.062 

82400 

6.01 

.63600 

1.103 

87308 

5.89 

.60129 

1.131 

91254 

5.874 

.58471 

1.456 
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A 

B 

C 

D 

X/D 

M 

Theta - 0 

Free Plume 
Expansion Radius 
(R/D) 

.96408' 

5.840 

.56404 

1.182 

1.05794 

5.80 

.53164 

1.242 

1.09118 

5.78 

.51986 

1.260 

1.14699 

5.75 

.50391 

1.292 

1.20885 

5.723 

.48549 

1.323 

1.27162 

5.715 

.46950 

1.355 

1.33487 

5.710 

.45436 

1.388 

1.43185 

5.69 

.43380 

1.437 

1.52983 

5.69 

.41471 

1.488 

1.68460 

5.69 

.38945 

1.561 

1.79133 

5.70 

.37429 

1.609 

1.90591 

5.705 

.36038 

1.660 

2.00453 

5.705 

.34923 

1.710 

2.10755 

5.715 

.33942 

1.743 

2.20516 

5.718 

.33107 

1.783 

2.31099 

5.718 

.32418 

1.825 

2.39077 

5.716 

.32047 

1.851 

2.40765 

5.701 

.32655 

1.862 

2.51795 

5.705 

.32794 

1.902 

2.6743 

5.717 

.31689 

1.955 

2.82225 

5.72 

.30810 

2.005 

2.96077 

5.722 

.30217 

2.050 

3.08055 

5.725 

.29927 ' 

2.086 

3.19207 

5.719 

.29961 

2.117 

3.35044 

5.685 

.31006 

2.165 
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It has been decided to break up the plume into increments at 
the following X/D locations: 

.56, .60, .7, .8, .9, 1.0, 1.1, 1.2, 1.3, 1.4, 1.5, 

1.6, 1.7, 1.8, 1.9, 2.0, 2.1, 2.2, 2.3, 2.4, 2.6, 

2.8, 3.0, 3.2, and 3.3. 


In addition, the graph shown in Figure 12 has been made and 
the values shown in Figure 13 have been chosen as the maximum Thetas 
and Mach numbers. With this data, the first part of the program may be 
loaded. 


1. Block 1 


a. There are 32 figures in each column of Table I. 
The value loaded here is X32. 

b. According to the X/D locations shown above, the 
region is to be incremented at 25 X/D's. Value 
loaded is X25 . 


c. Oblique shock theory is desired. X + 1. 

d. One half the distance between rocket engine center- 
lines divided by the exit diameter is .92773. 

e. Value is 43.004899. 

f. Value is 5760,0. 

g. Value is 1.23. 

2. Block 2 


The pressure ratio for this run was .0185. 

3. Block 3 

Load here column A of Table I in order going down, 
taking as many cards as needed. 

4. Block 4 


Load here column B 


5. Block 5 

Load here column C. 

6. Block 6 

Load here column D. 

7 . Block 7 

Load here the 25 values of X/D locations listed on page 28. 

8. Block 8 

Load here column A of Figure 13. 

9. Block 9 

Load here column ,B of Figure 13. 

10. Block 10 

a. The exit diameter of the engine is 6.666 feet. 

b. Value is .63. 

c. Since the last value in column A of the input 
table is 3.35044, the value loaded here is 
3.45044. 

11. Block 1 1 


a. Referring to Figure 16 and remembering to non- 
dimens ionalize by division by the exit diameter, 
the value loaded here is - .81404. 

b. From the same figure, the value (in radians) loaded 
is 1.57079. 

c. Value here is '2.09559. 

Now another table is made. Referring to Figure 17, 
locate the blockage regions (which are the heat shield 
and the five motors) and prepare the following table. 
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TABLE II 


Blockage 

x-coordinate 

A 

y- coordinate 
B 

z-coordinate 

C 

Radius 

D 

Heat Shield 

0.0 

0.0 

-3.4 

10.25 

Engine I 

-6.19 

6.19 

0.0 

3.33 

Engine II 

6.19 

6.19 

0.0 

3.33 

Engine III 

6.19 

-6.19 

0.0 

3.33 

Engine IV 

-6.19 

-6.19 

0.0 

3.33 

Engine V 

0.0 

0.0 

0.0 

3.33 


Also it is noted that the intersection region is located 
at X = -6.19, y = 0.0, Now this new data is loaded. 

12. Block 12 

a. There are six blockage regions so X6 is loaded here. 

b. -6.19 is the value as noted above. 

c. 0.0 is loaded here. 

13. Block 13 

Load here column A of Table II. 

14. Block 14 

Load here column B of Table II. 

15. Block 15 

Load here column C of Table II. 

16. Block 16 

Load here column D of Table II. 

17. Block 17 

a. Since a new intersection region is to be analyzed 

the first tendency would be to load 0.0 in this place. 



However, it is apparent that actually the new 
intersection region is formed by two motors separated 
by the same distance as before and formed by engines 
with the same characteristics as before. Thus, in 
effect the same intersection region is being analyzed 
but from a different point. 1.0 is input here. This 
tells the machine that it need not calculate new 
temperature and pressure profiles and that it may 
start execution with Chain II instead of Chain I. 

b. 1.0 is input here because clearly the intersection 

region is at a larger x-coordinate than is the point. 

The input for the first intersection region analysis 
is now complete. Directly behind it the data for the second intersection 
region is loaded. Of course, since the machine already has in memory the 
necessary intersection region data. Blocks 1 through 9 need not be input 
again. Starting with Block 10 and proceeding in the same manner as before 
the following is loaded: 

18. Block 10 

a. 6.6667 

b. .63 

c. 3.35044. 

19. Block 11 

a. - 1.97363 

b. .4895 

c. 2.09559. 

20. Block 12 

a. X6 

b. 0.0 

c. +6.19. 
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21 . 


Block 13 


Column A, Table II. 

22. Block 14 

Column B, Table II. 

23. Block 15 



Column C, 

Table II. 

« 

CM 

Block 16 



Column D, 

Table II. 

25. 

Block 17 



a. Load here now 1.0. The machine will continue to 
execute, will try to read data for Block 10 and 
finding none, will stop. 

b. 1.0. 

Figure 18 presents the input for this problem as it appears 
on the input cards. 


IX. CONCLUSIONS AND RECOMMENDATIONS 

1. This program may be used to estimate radiant heat flux rates 
to various missile bases from intersection regions produced by H 2 -0 2 
engines . 

2. Possible inaccuracies can exist in the calculated intersection 
region properties, the emissivity data, or in the application of equa- 
tion (15). 

3. Other inaccuracies arise because of the fact that second and 
higher order intersections are ignored, radiation from plume to plume 
is ignored, and blockage by intervening gas masses is ignored. 

4. Further work is needed in all the areas mentioned in 2 and 3 
above. In particular, an accurate method to calculate the region pro- 
perties and an accurate method of accounting for second order regions 
is essential before accurate estimates can be made. 



5. Measurements, preferably inflight, must be made before these 
analytical approaches can be completely trusted. Instrumentation must 
be developed that will measure these low transient heat flux rates. 
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FIGURE 2. GEOMETRY OF THE CROSS SECTION OF THE INTERSECTION REGION 
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OF THE INTERSECTION REGION CROSS SECTION 





Axial Line For Intersection Region. 
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FIGURE 5. SKETCH OF TYPICAL INCREMENT RADIATING TO A POINT ON THE BASE 





O GO 


TART , 

— — *— 5 L EM D A S H * 3 5011 Q 13 F O RTRAN 0 4 0 - 11 

1 COMMON L* GAMMA » XD ( 6Q ) *FM ( 60 ) * T ( 6Q i *RD ( OQ ) 

— C O MMON P-160) »TEnP 4 6U) »m)LS( .60) »AUIS(6Q) 

DlMtNSI ON A ( 6u ♦ 4 ) * T AbL ( 1 U ) » T Ab ( lu) * C ( 6U J 

3 RO) >LL »I. »MCUN »R.»PU »TU » GAMMA 

HCD »PBPX 


R CD» 4 ( A < I » J ) » I =1 » LL) » J = 1»4) 006- 

RCD»(XD(J)*J*1»L) ' ^ 007 

RCD-#-XA6J 006^ 

RCD » T Ab ' „ 009 

Qs^U , 010 

PRINT 1 ( Oil 

_L FORMAT ( 1H1 « 2 5X 1 7HRAD UT I ON PROGRAM//! QJ~2_ 

PRINT 2 

2_ FORMAT 127X9HCHA1N U)///) T 

•PRINT 1002 *PBPX 013 

1002 FORMAT ( 33.X-6HP6PXa t U,3».6/V-j 

DO 204 J= 1 t L 013 

AROgXD(O) 014- 

- DO 16 1 = 1 * LL ‘ 017 

' DAIA=A (.!-»- 1 ) Q 1 6 

IF (DATA-ARG) 16,4V»19 019 

16 CONT I NU E 1 Q20 

PRINT 216 » ARG 021 

GQ - XQ 204 022 

19 IF ( DA T A-A (1*1) ) 2U»2u.22 023 

20 PRINT 216 » ARG ' : T 1 : Q24 

GO TO 204 023 

22 D AXAEjlAX Xjb 1 . » 1 J : 1.Q 2 6 

F UNO { ARG-DAT AP ) / < DATA-DATAP ) 02 7 

EMHjslAX 1+Z-) — 02-6— 

FML=A(I-1#2) 029 

IF -4-EH HrFM L ) 3 0 1 30 > 2 7 : — - 03 0 

27 FMC=Fl)NC* (FMH-FML ) 031 

FM1 J)=FML+FMC , 1 062 

•, GO TO 32 036 

30 F44C=FUAICiU FML-^-fMH ) T Q34 

FM ( J ) =FML-FMC 036 

32 TF4=M-U 3 ) , -1— 0^6- 

TL=A(I-lf3) 037 

I F-t I T L4 38 3 6 3 3 066 - 

33 TC=FUNC* ( TH-T L ) N 039 

TU) = TL+TC - — — 1 T — Q4 0 

GO TO 40 041 

, 36 TOF U NC*(T L -TH) Q 4 2 

T ( J ) = T L-T C 043 

40 RQFUAt I »4.) : 044— 

RDL=A ( I -1 >4 ) 043 

llF (RDH-RDL ) 46»46»43 - r : 1 Q46 

43 RDC=FUNC*( KDH-RDL ) 047 

^,-RU-Uj.g KULdbKDC - ■ , . - - ,046 

GO TO 46 ' 049 

46 RDC rFU NCft I RDL.- R DH I 1 030- 

RD ( J ) =RDL-KDC 031 

46 GO IQ— 32 Q52 


FIGURE 8. RADIATION PROGRAM PRINTOUT 
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49 

pM( 1 ) = A f T - ? 1 . 

Q5-3 


T ( J 1 =AC 

054 


ij£LLJj_=Al T f 4 ) — 

0^6 

52 

I F ( 1.6-GAMMA) 61.63.63 

063 

6 1 

DR TNT 'J J r\ *(\A MMA * . 

064 


GO TO 71 

066 

63 

hn 7n ifsi.in ■ 

l 


IF< TABL(K.)-FM( J) > 70.7001.7001 

068 

70 

rnMT I Mi IF . . 

07X- 


PRINT 221 

072 


-GO TO 72 

07 3- 

7001 

IF ( TAB ( 1 C) -T ( J ) ) 71. 7b, 76 . ■ 

* * 074 

71 

DBIWT _ . 

Q..7-Q 


PRINT 224 

076 

.. 

tp/ i-l) 78*78*72 


78 

T ( J ) -A (1*3) 


72 

ip iMrnwi 74 . 74-73 _ 

077 

73 

MCON=G 

078 

7 4 

pp l T A = 1-1 H44 % 

0X9- 


GO TO 145 

080 

7A 

AA = (iAMMA+l • _ 

061 


AC=FM(J)/Q 

08^ 


AAA=OAMMA-l - 

QR X 


BB-FM ( J ) **U/G 

084 


00 = 1 t flAMMA#FM* 

0 ft 5 


D£LI=ASINFi l./FMIJ) ) 

• 086 


DLM=f f* ( AA*AC**L/-1 .+( AA*I 1 . +AA A*bti+A A*AC **4 * ))**-5 ) 

0 ft 7 


DELM=A$INF(t>LM**.5) 

088 


OL>tL = DELM-D£L I . 

n k w 


GO TO 100 

090 

99 


09 I 

100 

DELTA=DEL I +DDEL ' 

. 0V2 


THE T = T ANF I DFL TA)*t ( AA*riH) / ( ( FH { J ) *S I NF ( Dt L T A ) )#*G-1 -:)-!.) 

0 y 3 


. THE T A= AT ANF ( 1 . / THET ) 

094 


iFiArt^FiTHFTA-mji-^nnnn i45 1 iA5*i>n _ .: 

0 g 6- 

120 

GO TO 99 

096 

1 4 5 

IF (Mr ON) 145 f 157 *174 

Q9-7- 

146 

IF (J-i) 147,147*148 

098 

1 47 

PRINT ^07 

099 

148 

CP=2.*SINF(T( J) )**u 

100 



- --- lo 1 


0PC= < GAMMA-1. >/<J 

102 


op 0 = 1 .+OPO*FM**Q » 

___ 1 Q 6 


OPt-OPD** (GAMMA/ (VjAMmA— 1. ) ) 

104, 


SAV=OPP/OPt » 

106 


P ( J ) =P0*S AV 

106 


TPMPi n =Tn*5A\/** f / f,AiViMA- 1 - TVi-.Aj«MA 1 

107 


PRINT 208, XD( J) ,FFi<J) ,T(J) tPU) , TEMPI J) 



GO TO 3 90 - - 

110 

157 

IF ( J-l) 158,158,159 

111 

1 5ft 

pp I NT y ) 0 

112 

159 

CP-2.*MNF(T( J) )**U : 

113 


opp = ( r p* ( i.amma /u’) *FM*f y . ) + i . ...... 

1 1 4 


0PC= ( GAMMA-1. )/U . 

115 


OPl)=l.+OPC*FM»*U J : 

-116 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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bAV=OPP/OPE 


ALPHA* (GAMMA-1# i / ( GAMMA+1 • ) 


/ v. r MC I t t i i ) i 


VA=1.+ALPHA 


TEMP ( J ) =T0* ( (l.+BETA)*VA+ALPHA*<6ETA+tiETA**Q> ) / ( ( 1 . +VB ) *VA+BETA+V6 

9*BETA) 

PM til 20BfXO(J) tFM(J) »T( J) »P( J) #TEMP( J) 

r.f\ TA 1 OA 


174 IF IJ-l) 175, 175.170 

-1.75 ejaiM-f. Z12 i : : 

176 X=(FM#*U)*l ( SINF ( D6L.TA ) )**ui '' ! 

Y=U*GAMMA .. 1 * 

2=GAMMA-1. 


222 = 22**0 


XX=l.+(2/Q)*(FM**2. ) 

TyTf-T/Tl / V Y 


TEMPf J)=T2TC*TG ‘ 


P2PC=P2P1/ (XX** I GAMMA/ 2 ) ) 

p<jj S po*P 2PC 

PRINT 2oa,XG(J) ,FM(J) ,T(J) .PIJ) .TtMPU) 

- 1 9 0 C ( J ) =P EL TA ^ U J) 

IF (J-i) 1V2.1V2.1V4 

i a v uhki i \ - / t ii m c i r t itii-tt-r.iw i i — a i i . l i 1 


GO TO 195 

- 194 tsDlliJJ =HP 1 5 ( J- 1 ) + TANFI CI J- I ) ) * ( XU ( Jl -^ XUl J - l ) -) 

195 APIS! J)=SURTF(ROf J)**U~R**U> 

204 CONT-I NUE 

205 PRINT 206 

A CniYMflT (lull 


207 FORMAT (10X50HIS6nTkOPH. R6LA T 1 UN6 0660 For TcMP ANU P CALCU6AT.//) 

... 2UB Format I 23XSHXUI JJ-4.tr 15. a/ .cXoH Fm ( J ) = . 1 1 a . a , CXPHT-I J ) = .6 1 3«5 . 

V2X0HPIJ) =,tl3.5.2X6HTtwP( J ) = » t i 3 . 5 / / / ) 


■ f i i T UM iTW 


211 911HCALCULATI0N/// ) 


213 911HCALCULATI0N/// ) 


217 96HXu( J)=,ti6.B///) 


219 96HXUC J)=,hl6.b/// ) 


220196HGAMMA=,E16.6//) 

.... Z21 - ~£ X C t L 0 b -W U O U t GrtOC*. U 

222 FGKMAT (2bft06LiUUE bHOOC CANNOT bt UStO/ 

. -.2 Z^-9^BBB^-l- T CH-JUl--R AN ) C -.l- N t^ - H UG0N I U l/ J - A l 

224 FORMAT U6HD£LTA = 6i> GtGREtb///) 
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FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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246 


309 ARt AC=0 • 

-344) AK fcA C * ARE AC + ( TEMP ( K . W ) ** 4 4* ( X -D ( (W ) -XD ( *4 

TOT A=ARt AC+CAREA 

— A VT £ M( J )=U ) TA/ 4 XDS=X0( J4 ) 

AVTEMt J)=AVTEM< J)**.25 

r.r\ m -a i 


314 M=l-1 

DO 31.6 XsJ.M 1 — 

IF ( K.-J ) 313,317.310 

317 AKtACaq. . 

318 AREAC=AREAC+(TEMPIX)**4)*IXD(X+1 ) — XD t «. ) ) 

T O T AsA R EA C 

AVTEMI J)=T0TA/ (XDI I )-XD< J ) ) 

A VT t Ml J )» A VT t MU ) *IU 25.. 

321 IF IXDIU-XDS) 322,331,331 

-32-2 CPKESsIPtL) ) * tXUO-X-01 L ) -) 

N = L-1 . 

nn 1^7 ,u> l.u 


IF IXW-J) 327,326,327 

326 — A R EA Pg Q. 

327 ARtAP=AREAP+IP(XW) ) * ( XD ( XW+1 ) -XD ( XW ) ) 

t Q IP = ARE A P+ C P R E6 

AVPI J)=T0TP/(XD8-XD(J) ) 

60... T O— 330 

331 MN= I -1 

PI5-335- X = J-»W 

IF IX-J) 335,334,335 

_334— ARSAP^q, 

335 AREAP = AREAP+(P(M ) * I XD I X+l ) -XU l X.) ) 

T Q TP = A RE A P - 

AVP(J)=TOTP/(XD(I )-XD(J) ) 

-558--- P E ED 1 J ) s -t-F-L E I J ) ) *D 1 A 

6I<j2=1./SURTF(FX**u+( 6XD + XD(J) ) **u ) 

5T-05. S l. / SU RTF I FX**U + I 5 XD+ XD I J + l ) ) »» 0) 

XXX=SIG2*(-U*ATANFITDIS< J)*51G2) ) 


F(J)=(FX/6.283185)*<XXX+YYY) 

PH 2U (J) = FX I » A VP I J ) 

PL=PH20(J)*OLED(J) 

1- FIP L - .Q 0 1) 340 7,3 3 5 0 ,335-0 

3390 IFIPL-.02) 3391 » :1394 . 3 394 

.3391 IF ( AVTEMI . I >-2200. ) -3592 ,-3392 ,3393 — 

3 392 Ap = .691BE-2-.10 35E-4*AVTEM< J) + .5906E-8*AVTEM< J)**2 

A1=4.1462-.3099E-2*AVTEM( J )+.9596t-6*AVTEM( J)**2 

... 9-_>l 2 1 E--3»A V T E M. ( J ) » * 3 - 

A2 =-4 1. 608+. 2 765E-1* AVTEMI J ) -.243t-5*AVTc-Ml J)**2 


A3=175.61-.8474E-1*AVTEM( J)-.2372t-4*AVTtM( J)**2 
-9+-. 1 303 E-7*A-VTEM I J ) **3 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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EMISI J.)J=AO+AlttPl.+A?*Pl ,»*?+A3«Pl «*3 

GO TO 3406 

3393 A0=.1042£-2-.95B2E-b*AVTEH( J ) +. 1 574E-9*AVTEM I J>»*? 
A1=2.524-.10BBE-2*AVT£M< J)+.1197E-6*AVT£M< J)**2 

A7=-7H.9A +. 1 19 4 A t- - 6* A VT i- w ( ,j 1 *tt P 

EMI S I J ) =A0+A1*PL+A2*PL**2 

GO TO 3406- 

3394 IFIPL-. 1)3395. 3395. 3400 

33.95 I F ( AVTEMI J 1— 2200. ) 3396.3396.3397 

3396 GO TO 3392 

3 19 7 I F ( A V T t Ml— LI -.340 tU- L -3 3 9G . 3 3 9 9 . 3 1 9 9. 

339B A0=.569E-2-.2627t-P*AVTtM( J ) -.1168t-9*AVTEM < J)**2 

. 9+.B381 E— 1 3*AVTEM( J ) **3 1 

A1=2.594-.1468E-2*AVTEM< J )+.3207E-6*AVTEM( J )**2 

9- . 2 583E — 1 0* A VT E M I J ) ** 3 

A2=-.7B73-.2297E-2*AVTEM( J)+.923BE-6*AVTEm( J)**2 

9- . 76 0 5 E -10* A VT E M (. l ) **3 

.A3*-.8534E+2+.845Bt-l*AVTEM( J ) -. 2441E-4*AVT EM ( J)*»2 

9+«7 1 19E.-B*AVTEM I J )**3 , _ 

EM1M J)=A0+A1*PL+A2*PL**2+A3*PL**3 
GO TO 3-406 : 1 Lj 

3399 A0=.1052-.731BE-4*AVTEMl J ) +. 1645E-7*AVTEM I J)**2 

- 9-.170HF-1 l«AWTFMt.n»*3 i 

Al=-2.299+.2014E-2*AVTEM< J)-. 4971 E-6*AVT EMI J)**2 

9+. 3741 E-10* A VTEM ( J ) ** 3 

A2=66.19-.4732E-1*AVTEM( J)+.1086E-4*AVTEM< J)**2 

9-.79B6E-9»AVTEMU)»»3 _____ 

A3=— 405.9+»2897#AVT EM I J ) T.67ult— 4*AVTtMI J ) **2 

9 + .5014£-B»AVTEMI.J)»»3 ■ 

EMI 5 1 J)=A0+Al*PL+A2*PL**2+A3*Pl_**3 
: GO TO 3406 _ 

3400 I F I PL-1.0 ) 3401-.340l.3402 

3401 A0=.1Q06— .1014E-3*A\^rFM( J 1 + .439QE — 7 + AVTEM ( J)**? ; 

,9-.9089E-ll*AVTEM{ J ) **3+. 709 bE-15*AVTEM < J)**4 

Al=.77Q5+.lnl6F-3*AOTFMI J ) -. 7/HBF-b»AVTFM I J ) »♦? 

9+.6294E— 10*AVTEMl 3 ) **3~. 5309t-14*AVTEM I J ) **4 

_ A2 =-1.17B4-.3062E-3»AVTEMl J)+.5419£-6»AVTLMI J)»»2 

• 9-.1516E-9*AVTEM< J > **3+. 12B9E-13*AVTEm I J ) **4 

A3=.'9306+.P469E-3*AVTEMI J ) — . BbB 5E-6 + AVT EM I J)*»? 

9+» 1773E-9*AVTEM ( J ) **3-. 1454E-13*AVTEM ( J ) **4 

A& = -.77 17-.2q7qF-3*A\/TFMl .))+.7Q3H E-A»AU7H» l l 

9-.7382E-10*AVTEM< J)*»3+.5799E-14*AVTEM(J)**4 

EM16I J I=AQ+A1* P L+A2* P I. ** 2+A3*PL **3+A4* P L** 4 

GO TO 3406 

3402 IF IPL-IO. ) . 3403.34^3. 3404 — — — , — 

3403 A0=.28l4-.t>083t-4*AVTEMI J)-. 99BE-8* AVTEMI J)**2 

9 + .26B4E — 1 1«AVTFMI.))»»3 , 

A1=+.9226E-1-.3216E-5*AVTEM(J)+.1709E-8*AVTEM( J)**2 

p 9-.B164E-1?*AVTFM(-J)«»3 1 _ 

A2=-.1176£-l+.3U54E-5*AVTEMl J)-.15i4E-8*AVTEM( J)**2 

9 + .2766E-17»AVTFMl J)»»3 

A3= .5398E-3-. 1976t-6*AVf EM ( J )+ . 9794E-1 0*AVTEM I J ) **2 

9-. 164AK -1 4 

EMI6I J)=A0+A1*PL+A2*PL**2+A'3*PL**3 
GO TO 3406 ! T — 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued! 
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IF (J-l) 350*3 46,330 

_ 348 — QRADS = Q~» _ — 

3SF=0* 

— 3EQ — QRAD3=SAVE+QRADS 

SSF=F(J)+SSF 


QR AD ( J ) = 3 A V£ 

—PRINT- 4 27 *-QHA D, ( J ) »URA DS : 

42 3 FORMAT ( 3HdLtD= * t.13 • 5 , 2X6HAVTtM= , L13 . 3 »2X3nPH2Q= »L13*5* 

. 9 2 X4 HSS Fi= ,EI3. 5 » 2 X 3n E M iS = »bE3.« E// ) * — 

42 7 FORMAT ( 15X3HURAD= * Li 3 . 3 , 4X6NuKAD3= » tl3 .3/ ) 

-A 40 CO NTIN UE < 

PRINT 426 

4 26 F O RM A T i 1 H -*h 2 0X 21 HU b Y AVT - + AVP METH OD //) — 

PRINT 408 3 *GRAD3 

40 83 FUK MAX-U.7iiXUlAl — KAD l rtNT H b AT F lU X 16 1 1 6 • o » 2 X 1 3H 6 T U / 3EC - 3W FT//) 

414 PRINT 415 

,4,13 FORMAT < 1H1 ) 

430 FORMAT (10X50HPL 13 LL3S [HAH #001 FT-ATM - PROGRAM WILL DISCARD//) 

43 1 F GKMAT (lOX 49 HPL GRlhILR T-HAH- 1-QO-F J-^ATM - PROGRAM WILL DI8CAKU//0 — 

CALL CHAIN (3) 

— END - - - - ■ ■ - 

* FORTRAN 


COMMON P ( 60 ) *TEmP( 6(J) * dD I 3 l 60 ) » AD I 3 t 60 ) 
COMMON FT 6.0 J-*- ANGLE L6U ) »FLt(6U) ,P-HLul6U-l»bm: 
COMMON FR»3XD,ORAD( 60) »OKL>A(60) 

.-DIMENSION MA H 51 ) — 

DO 416 J= 1 * L 


351 PRINT 332 
PRINT 4082 

330 ... ABSi 3= (FLt< J) i/30. 

CN=Ab3 I S*3 1 NF (ANGLl( j) ) 

D Q i9 60 1 =1, 51 

IF (1-1) 361*361 *365 

341 CR=XDL3) 

C I = I ' 

3 A d 3 L3= Ad414* C 4 *01 m- 

GO TO 369 

-363 C 1 = 1 

AI = U-i. 

3 Ad 3 1 3=4631 3-* Cl *0lA 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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fc. 


CR s CN* A l+XD (J ) . 3Q7 

369 DO 392 K= 1 ,L 306 

- — lEICR-XOIKI) 374.371.392 : 3 p 9 

371 TW=TEMP ( X ) 310 

-3XX- 

GO TO 393 312 

3 74 FUNT S I CR-XU I K - 1 ) ) L I AbSF I XDI M- X U U-U ) ) a | a 

QTH=TEMP l X ) 314 

_ QJL = TEMP(K.-1 ) • 3X3 

IF (UTH-QTL ) 381.381,378 316 

338 TF = F|,INT* (.uTH— UTL ) 3 1 7 

TW=uTL+TE 316 

GO TO 383 * ; 3X3— 

381 TE=FUNT* I QTL-UTH) 320 

TW=OTL-TE 321 

383 QPH=P(K) 322 

— : QPI =P1 K-l ) 3X3— 

IF(OPH-QPL) 389.389,386 324 

336 pc= Flint* 1 oph— opl i 329 - 

PW=uPL+PC 326 

GO TO 395. 327 

389 PC=FUNT*( UPL— OPH ) 326 

PW = OPI -PC ! : : 324 

GO TO 395 • 330 

-3 32 CONTINUE ?? j 

TW= TEMP ( L ) 332 

P W=P 1 L ) 1 333 

395 POW=FXI*PW 

POWI =POW*.5Ari.S 1 ,s 

3950 IFIPOWL-.02) 3933.3954,3984 

3953 IF ( Tw-2200. ) 3956.3956*3983 1 

3956 A1=4.1462-.3099E-2*T*+.9596E-6*TW**2-.121E-9*TW**3 

A2 =— 4 1 .808+.2 7 65F— l* T ri-»2 43L-5*TW **2-»10t- B * T W ** 3 , 

A3=175.61-.'8474E-l*Tw-.2372E-4*TW**2+.1303E-7*TW**3 

VALI I )=A1 *PQWt2.»A2*IPQd**2 1 *.3AtL5 1 6 + 3 . *A 3* I PUw**3 ) * I aABSi 6**2 ) 

GO TO 3960 

3957 A l = 2.524-.10B8E-2«T> l H-.119 7b-b»Tw**^ 

A2 = -28.96 + .1552£-1*Tia/-.1945E-3*TW**2 

, VALI I l=A1*POUH-2.»A2»(PuO*»2)».SAHSI.S 

GO TO 3960 

3954 I F ( POWI -.11 3959.3959. -1961 1 

3959 IF I TW-2200.) 3962 » 3962 » 3963 

3-962 GO TO 3956 ! 

3963 IF I TW-3400. ) 3964 » 3 963 » 3903 

39 6 4 A I = 2 . 3 94- . 1 4 6BE - Z *TW+ . 3l Q 7 E ~6* T W* - * 2 ~ . z 5 8- 3 1 - 1Q *T W ** 3 I — , 

A2=-.7o73-.2297£— 2*T*+.9236t— b*TW**2-.7603E-10*Tvl**3 

A3 = -.8 534F + 2 + .a45Bl--l*T»i-.Z44 I F -4»T w** 2+ . 2 I 1 QF-H*TW**3 

VALI 1 )=Al*POW+2.*A2*tPUW**2)*6A6i>I8+3.*A3*< POW**3 ) * I 8A6SI S**2 ) 

G O T O 3960 

3965 Al=— 2 .299+.2014E— 2*TW— .49 71t-6*TW**2+.3741E-10*TW**3 s 

A2=66«19-.4732£-l*Twt»-Lfl66t -4*TW* . * 2-.79832- 9*T W** 3 „ 

A3=-403.9+.2897*TW-.e70it-4*TW**2+.3014£-o*Tw**3 
VALI 1 )=A)*PlJU + 2.*rt2*IPw.l*»2)*.SArtM.S+3.»A'> »IPij|»»»3>»I.SAn.sl S**2) , 

GO TO 3960 

3961 IFIPOUiL-l . D) 3966.3966. 39a7 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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£-966 41 ^ -7 E -=-3*.T Vj ^ ^£-£6 bfc^ fa*Tw**2+ . fe294t-lQ*Tw** £ — 

9-.6309E-14*TW**4 

42^-- 1^447£ 4^-30 £2E ^£*4^- ♦-£6-19 t - r lsit>E-9*Tfr**3 

9+.1289t-13*Tw**4 

A^=-^ 3 ilh+- t i 4 K 0 i--^*Tui-..ARn»-ii--ri*Tl.i*»^ + T t 7 7' 1 r-V*Tft**3 — - 

9-. 1464t-l3*rw**4 

— A 4=- « 2 717- *24Oi^.£iUU6*»2 9 bttt - fr M'W**2- » 7 £tt2fc - l 0*T w **3 

9+.b799E-14*TW**4 

- V4-U I 4— A-l *P 0 W-+ 2-*JfcA2 *4-P Q#**Z4Jk± A #£ l - b+3 ♦ JM4*W*- *3 ) * ( S ADb4b **2 ) 

9+4.*A4*( PQW**4) *(bAobIS**3 i 

GQ T O 3960 

3967 IF(P0WL“10*) 3 966 > 3 966 * 396 9 

£9 66 4U- ? + # 9 22 6 6 -1^3^64-d* Tav * ♦! 7 Q 9 L-^ 6 *-Uv *^^*6 1646-1 2*T 

A 2 = “.il 76 E-l+# 3034 t- 3 *Trt-.i 3 l i +t“ 0 *Trt**<' + #i '7 66 t-it*Trt**:> 

A 6 = . 3.3 £ tt E -^-.«£ 9Z6^-6 *-T T£*-*2-^ Xo 4 6 t - 1 6 *4 ^**6 

V AL ( I )=Al*P0w+2.*A2*tPuw**2)*6Abbl6+3.;*A3*{P0***3)*(bAbbIb**^ ) 


3969 

CQ TO 396 Q 

Al = .6669t-2-.bl:>lt:-&*Tw-.9l30E-lO*Tw**2 


4Q6 0 

V AL ( I )=Al*P0W+2.*A2*(P09\/**2 )*bAcbl6 

VA1 ( I ) -VAL ( IT *TW**4 



DO 400 M = 2 > 60 > 2 

336 


IL (M- 2 ) 'avjw f C 190 r 4 ( p M» 

337 

399 

W E = 0 • 

338 

400 

WF=Wt+VAL(M) 

339 


DO 4u4 M=1 >61 >2 

340 


IF(M-l) 403 >403 > ^04 

34 j 

403 

YOU= 0 . 

342 

604 

YOU= YOU^-y A1 (MJ 

44o- 


F = F ( J ) 

bA V 6 = Ad bib* 4.7611 6 " 1d*F* 1-UAL U 4-+6* **t-+2 . *-Yuu-VAL ( 6X41/ 3* - -34A- 

IF(J-I) 406 >40 7 >406 346 


4U7 UKUAS=0* • - — 34b 

4U8 URDAb = bA*V t + UkUAb 

PR, Ik T 426>J»X0UM 

426 FUKMAT ( 3UX2hJ= > 13/ >22 XoHaD( J ) = ftl3«6// ) 

URlAI J ) =bAVL ... — - 


START 

CHAIN 

4 


PRINT 427 > UKDA ( J ) > OKDAb 

627- . FORMA Li lbX6FiUR4^ »LU. t > »4 X6H u KU A b = ♦ El - 2 > 3/ ) — — 

416 CONTINUE 

PRINT 6062 . 

4082 FORMAT (20X30HHEAT FLUX bY AobUKPTlON Pit THOU / ) 

MUN-I- ^40££ » URO Ab — — 

4083 FORMAT ( 27HTOTAL KauIANT H t a 1 FLUX lb cl6 .6 , 2Xl3HDTU/6LL-bO 

- .-PRINT Alb - — 1-4 

416 FORMAT ( lhl ) 

CAUL CHAIN (4) : 

421, CONTINUE 

• • - END— - — - ..... — - - - -• - • - - 

* FORTRAN 

— DIMENSION- Ai_40-»44»TU1£(601 - — - - 

> COMMON L >0AMMA >XD ( 6u ) >pM ( 60 ) * T ( 60 ) #Ko ( 60 ) 

CUtuYiUN , P L 601 » J-LmR Loo ) » oU lb LbU4 » ad Ibl 60 X * 

COMMON ' F ( 60 J > AiMOLt ( 6U ) , FLt ( C>0 ) >Pn2U < 60 ) * DMA 6L tL) 1 A >Pn i > FX1 
COMMON. -Fk »4>XU+QKAU ( 601- - » utkUAlbQ 1 _ ._ 


FT// ) 


£4£ 

349 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 
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, — Sf.U . 1. 1. L L ,..F.l NtX»f.lft.tS 

RCD,(<AII»J>,I=1»LLLL),J=1,4) 

RCO.FLAG.PX ; 

. PRINT 252 

print thA : 

252 FORMAT (27X9HCHA1N (4)//) 

253 FORMAT < 1 7X27HriL0CXAGE CORRECTION FACTORS//) 
IF (PX) 32.32,100 

..100 SI-ASINFIAe)SF(FINTY/(FX»I1IA) )) 

XP= (COOK ( SI ) )*(ABSF(FK> )*0IA 

IPI XP-F IUTX1 30 . 3 1 .31 '. ' .. 

31 PX= ( XP-F I NTX ) * ( _ 1 • ) 

GO T-Q --32 

30 PX=IXP+FINTX)*(-1«) 

-32 P-Y-=o« : 

P2= (-1. )*SXD*DIA i 

110 69 Ls-L+J 

Q = 2. 

V.G g -SlN FT PHI ) * »Q ». : ; 

VGG g COSF (PHI ) **U • ' • 

VC- AO 13(J)*B01S. 1J) , . 

VO g ADI S ( J ) **Q 

VE=HOt.S.tJ.U*.Q 

TDIS g VC*ISQRTF(l./(VE*VGG+VD*VG) ) ) 

RC=VC* ( SURTF < 1 «./ ( Vt *VG+VD*VG0 ) ) ) : 

XCOR g TOIS*(SINF(SI > )*0IA 

Y£OR=IO.IS*1CQG£.1G-I1.I»01A 

PRINT 999, J 

. 39,9,. FORMAT 1.2 AX/H.ls t.l 3// > ■! 

YR=F I NT Y-YCOR 

YL g F I NT Y+YCOR ■ 

IF(FINTY)26,27,27 

_26 XR g F I NTX— XCQR ’ : 

• xl=fintx+xcor 

fid TO 2ft : 

27 XR=F I NTX+XCOR 

XL =F I NTX— XCQR ... 

2d CONTINUE 

DP 4 .9. JC = 1 , i „3 ! - 

IF ( K.-1 >21,21,22 

2? I F ( X — 5 1 '2 1 ♦ / W4 v 

24 IF ( X-3 >21,25,25 

—21 PRINT-. 2. 9 . '.. : 

29 FORMAT (19HLtFT Slut OF REGION//) 

-99.8 FL.Y g O .O ! 1 : 3 

AINIX=XL-PX 

.. A 1 N.T, Y s.Y.I 

GO TO 40 

.23: PRINT .41 . 

41 FORMAT ( 16HMIDDLE OF REGION//) 

aintx = f:intx- p x ■■ ___ __ 

ainty=finty , : 

GO to 40 ; £ 

25 PRINT 42 

—42— FORMAT ( / OHR IGH.T— SIDE O F- REG I ON/-/-) 


FIGURE 8. RADIATION PROGRAM PRINTOUT (Continued) 






IF ( ) 70*71 >70 



7 LL 

J PFF=( s(j^TF LFiUt#; ^ ) ) 


HYP = ZO*<FFF/ (SXD+XDU)) ) 

At) TO 'Ll _ .. : 


71 

HYP=ZO*l ( Ab6F ( F tC ) -KC ) / ( iXD+Xl) ( J ) ) ) 



2uL 

XQf 1C0.S>F( THtTA) )*HYP 



Y0= ( 6 I NF ( THtTA ) )*hYP 

— i£— LAXN.IX1 30* 5 0^3 L ... . - 

50 XO=-1**XO ‘ ' ■' 

37 IF CAINTYi- 10 %lu il-i- 1 _ : 1 

10 Y0=-1.*Y0 

11 H = A t LI » 1 ) 

PRINT 10 02 » X0 * Y0 

1 UU 2— FORMAT t 3HXO- »E I6»fl * 2X3HYU= » 1 16 « b / / ) 

FFtC=A(LI*2> 

RAUlUS=A(Ll >4) — : 

EUN=(XU+PX“H)**2.+( YO-FFA)**2. ' 

SA 0 IUE sR A D IU6 ** 2. 

IF (SADIUS-EQN) 79.49,49 

- —79 CO NTIN U E — * 

1 FLY=1.+FLY 

AS CO NT I N U E _ — 

QRAU( J)=QRAD( J) *FLY/3. 

QRUA i J ) gQROA ( J ) *FL Y / 3. 

5 uRAU6 s uRAUi> + URAulJ) 

: OHOASg-URDAS^ORO^aJJ _ ■ - 

PRINT 1005,URAD(J)*uRADS 

PRINT 1UQ6 »ORDA-ULUu KPA S — - - 

1005 FUKi v iAT ( 5HUKAD= *E16.b ,2X6HuRAu>b= ,E16.b// ) 

10 06 Furm A T ( . frN u KO A ^ Llo^tt.^Xfc H u KU M^s *e.l6.a // ) 

69 CONTINUE 

4036 FORMAT (29HwITh di-UCNAbt: cuKkLLT i oN * Tht 

• 9A6HRA0-1 A-T4A/ L E - LU X- o Y THE AVT— - * NU * ¥F ♦ MfcTttOP IS Elto .b . 2X3 - H U- I - U 

910H/SEC-SU FT/ 

9 -32 H 0 Y THE AttSOB g - T I ON M E TH O D # I T 1 6 fcl fe. fl/// ) r 

PRINT 5000 

6000 FORM A T (lHl) — — - 

IF(FLAO) 419 #420 >419 

420 CALL CHAIN (4 -) — _____ ■ _____ 

00 TO 421 

449 CAL L CH A IN (2 ) 

421 CONTINUE 

END - - - - - — — — — 


FIGURE 8. 


RADIATION PROGRAM PRINTOUT (Concluded) 
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CHAIN 2 




►DO NOT SATISFIED 


-DO SATISFIED 


FLE(J)-25^_ 
so Lve ANGLE f .) 


FLE(J) *50 
ANGLE ( J) =90° 



FLE(J) “50 

i >0 



v l > 0 SOLVE 

~\ 




| "AjjirTJ 

; f * FLE(J 

lJ 



DO NOT SATISFIED - 


<JJD(I)-Xpj 

CONTINUE 



DO SATISFIED - 


SOLVE 

. 


sol vi : 

AREAC 


TOTA 


SOLVE 

AVTEM(J) 



FIGURE 9. FLOW DIAGRAM (Continued) 
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SAT 

ARSAP-0 

ARSAP-0 

I 



SOLVE * 1 


1 SOLVE 

A REAP 


A REAP 



t— 

SOLVE 


SOLVE 

TO TP 


:otp 

* 


J 

1 snt.vrc 1 

.J SOT. VS! L 

1 SOLVE j 

1 MZU'l — 1 ' Ui-iiliLU 1 l A7p(J) 1 


[ S0LV1 

C3 

Ixxx j 

[ 

SOI vv 

— tOJ 



|0j SOLVE 
H EMIS(J) 


SOLVE 

1 -Q 

^^lOoNv. 

>n 

' .1 

STOP 

EMIS(J) 


I SOLVE | 

L 

1 

ERR. 


DO NOT SATISFIED 


> *° r 

i/v nc«o 

^ 

8SF-0 


FIGURE 9. FLOW DIAGRAM (Continued) 
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FIGURE 9. FLOW DIAGRAM (Continued) 










m 





YOU-0 


gUADS-0 


FIGURE 9. FLOW DIAGRAM (Continued) 




























FIGURE 9. FLOW DIAGRAM (Continued) 


























































FIGURE 9. FLOW DIAGRAM (Continued) 







RADIATION PROGRAM 
CHAIN tit 


PBPX* 0.185OOE-O1 

OBLIOUF >HOC-< RELATION U5gh FOR TEMP AND PCESSuRE 

calculation 





XOIJ1* 

0.560006 00 






PrtI Jl 


0.65*7 6 

01 TIJ) 

* 0.600R6P 00 

PU| 

0 . 204216. !I0 

TEMPI J)» 

0.331616 

04 




X D | J J ■ 

0.6P0O0E 00 






FMt Jl 


(1 ,S6N*' 6 

01 TIJ) 

* 0.61031P 00 

PC J| 

0 .2092QE- 0 

TEMPI J|* 

0.335696 

04 




X D 1 JI* 

0.700006 00 






PHI Jl 


0.M041E 

01 TIJ) 

* 0 . 69 0 9 4 P no 

PIJI 

0.249936-00 

TEMPI JI* 

0.377206 

04 




XD[ J). 

0 « 0 0 0 0 OF 00 






PH1JI 


O,*0l5<*6 

in Tui 

« 0.64358F 00 

PC JI 

0 .222996- io • 

TEMP (Ji* 

0.350756 

04 




XD | J ] ■ 

0.90000E 00 


0.195276-00 




PHI Jt 


0.5P791E 

01 Ttjl 

■ 0 . 5699 8 c nO 

PIJI 

TEMPI j|* 

0.32248E 

04 




X D | J I * 

0.10000F 01 






PHIJJ 


0.5*?476 

01 TIJ) 

■ 0.55164F no 

PC JI 

0.175976-00 

TEHPf JI* 

0 . 31)2786 

04 




XD| J| • 

0.110006 01 






FMCJ1 


0.5775 'E 

01 TIJ! 

• 0.517T4F 00 

PC JI 

0.159436-00 

TEMP | J I ■ 

0.2*5*76 

04 




XD ( J 1 * 

0,120006 01 






PHI Jl 


0.57?6<-E 

01 TIJ] 

■ 0 . 4a8txc.no 

PIJI 

0.145976- ;0 

TEMPI JI* 

0.27210E 

04 




XD1 JI« 

0 « 130 0 OF oi 






pm I J 1 


0,5712*6 

01 TIJ) 

* O.4627iP-n0 

pc JI 

0.134456-00 

TEMPI j|* 

0.260306 

04 




X D [ J ) ■ 

0.140O0F 01 






FMt Jl 


0 ,5696*E 

01 TIJ] 

* 0.44085F-00 

PIJI 

0 .124016-tjO 

TEMPI J|* 

0.250416 

04 




XO( J!« 

0*150006 01 






F M 1 J 1 


0,56R0'iE 

01 TIJ) 

■ 0.42059P-10 

PC J) 

0 .116326-uO 

TEMPI J|* 

0 ■ 241686 

04 




XO| J|- 

o.monOF oi 






PMtJI 


0 ,H6«on6 

Ot TIJ) 

* 0.40396F-00 

PIJI 

0.109196-00 

TFMPf JI* 

0.234356 

04 




xn ( ji * 

0.17000F 01 






*P M t J 1 


O.H69l-tE 

01 TUI 

■ 0.367?6F«n0 

PIJI 

0.10277E-O0 

TEMPI Jl* 

0.22773E 

04 




X 1) f J I * 

0,100006 01 






PHI JI 


0 , 570 0 •*£ 

01 TIJ) 

* 0 . 373?4e-00 

PC J) 

0.972416*01 

TEMPI J)- 

0.222026 

04 




XD | J|« 

o.ioonOF 01 






FMt Jt 


0.S7047E 

01 TIJ1 

* o. 36iinF.no 

PC JI 

0 .9261QE-01 

TEMPI Jl * 

0.217236 

04 




XD | JI* 

0. 2000 06 01. 






PHI J| 


0,5705 E 

01 TIJ] 

■ 0.34974F-00 

PIJI 

0. *84196-01 

TEMPI JI* 

0.212*96 

04 




XO[ J] « 

0.21000F 01 






PMtJI 


0.8714 .16 

01 TIJ1 

a 0.34014F.00 

PC J) 

0 . *48806- Ul 

TEMP IJ 1 * 

0.209226 

04 




XD| J|« 

0 , 220(1 Of 01 






PHI J! 


0.8717-6 

01 TIJ] 

* 0. 33181F-00 

PC J] 

0 , *1 9066- >1 

TEMPI J)« 

0.206026 

04 

PHI Jf 



XD I J 1 * 

0.23000E 01 







0,8716 6 

01 TIJI 

■ 0.32A90F-O0 

PIJI 

0.795136-01 

TEMPI Jl* 

0.203646 

04 




X D | J | ■ 

0. 240006 01 






PHI JI 


0.87n7*6 

01 TIJ) 

■ 0. 32379F.nO 

PIJI 

0.792216-n 

TEHPf Jl* 

0.203336 

04 




xnt j|. 

0,2*0006 01 






PMtJI 


0.8711^6 

01 TIJI 

■ 0.32214F-00 

PC J) 

0.78626E--1 

TEMPI Jl* 

0.202716 

04 




XD I J ) ■ 

0 . 2*00 OF 01 






FMtJI 


0.871986 

01 TIJI 

■ o. 30942F.no 

PIJI 

0 . 742826-0 1 

TEMPfj}- 

0.198186 

04 

PH | Jt 



XD I J | ■ 

0.30000E 01 







0,8793 6 

01 TIJI 

• 0.3D122F.OO 

PIJI 

0.715706-01 

TEMPI J). 

0.195346 

04 




XDt J|* 

0 , 3?0 0 OF n 






PHI J! 


0.871 7 IE 

01. TIJ) 

« 0. 30013F.nO 

PIJI 

0. 712666-01 

TEMPI Jl* 

0.195026 

04 




XD| J|* 

0 . 330 0 OF 01 






PHI JI 


0.86958-g 

01 TIJI 

• 0 . 3O673F.O0 

PIJI 

■ 0.735936-01 

TEMPfJI* 

0.197466 

04 
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chain t2i 


0 BY AVT • A VP METHOD 

FK**0.81404f on 
PH I • 0 . 15708F 03L 

—Jm~ 1 " 

XD t J 1 « 0.86none 00 


BLED* 

0.?66?9F 02 AVTEM* 0.247636 04 PM2f)» 0.661076-01 

SSF* 

0.984206*06 

6MIS* 

0.228776-00 


QRAn* 0.40307E-05 QRADS* 0.403076-05 






J* 2 

xDui* o.^onooe oo 





BLED* 

9.87655E n2 AVTEM* 0.223446 04 PW2n* 0.549856*01 

SSF* 

0 . 16411E.Q4 

EMIS* 

0. 311776-00 

--- 

QRAO* 0.570746*04 QRADS* 0.61105P-04 

j* . 3 

XDtJJ* 0.700006 00 





BLED* 

0 . 1 2166E 03 AVTEM* 0.209166 04 PH20* 0.49935E-01 

SSF* 

0.60190E-04 

EMIS* 

0.396076*00 



QRAn* 0.158016-03 QRADS* 0.21912F-03 

J* 4 ' - ' 

X D ( J ) * 0,800006 00 


■ 

-- - - 

- 

BLED* 

0*1 9894F 03 AVTEM* 0.202966 04 PW2n* 0.481656*01 

SSF* 

0,129446.03 

EMIS* 

0.44783E-00 


QRAO. 0.250526-03 QRADS* 0.46964P-03 

J* 5 

XD [ J ] ■ 0, 900006 00 





BLFO* 

0.J3330E n3 AVTEM* 0.199946 04 pM2o* 0.47?486-01 

SSF* 

0.217486-03 

EMIS* 

0.484616-00 


QUAD* 0.324666-03 QRADS* 0.794306-03 






J* 6 

XDtJJ* 0.100006 01 





bled* 

0,333306 03 AVTEM* 0.19935E 04 PH20* 0.470936*01 

SSF* 

0.319326*03 

EMIS* 

0.484976-00 


QRAO* 0,371376-03 QRADS* 0.11657F-02 






J* 7 

X D ( J J * 0,110006 01 





BLED* 

0.166656 03 AVTEM* 0.20035E 04 PW2n* 0.475666*01 

SSF* 

0 . 431376-03 

EMIS* 

0.42864E-Q0 


QRAO* 0,368436-03 QRADS* 0.153416-02 






J* 8 

XD 1 J 1 ■ 0.12O006 01 





BLED* 

0.166656 03 AVTEM* 0.199706 04 PH20* 0.473506*01 

SSF* 

n, 550946. 03 

EMIS* 

0.42885E-00 

- - — 

QwAu* 0.388346-03 QRADS* C.19224P-Q2 

J* 9 ' - 

xDtjj* 0.13*006 oi 





— ■ - - -- 

BLED* 

0.166656 03 AVTEM* 0.19921E 04 PH20* 0. 471686101 

SSF* 

0.676026.03 

EMIS* 

0.426996-00 


QRAO* 0.402296*03 QRADS* 0.23247F-02 






J* 10 

XDtJJ* 0.140006 01 





BLED* 

0.166656 03 AVTEM* 0.19882E 04 PH20* O.470l5ElQl 

SSF* 

0.804966-03 

EMIS* 

0,429076-00 


QRAt)* 0.411566-03 ~ QRADS* 0.273636-02 
J* 11 

XDtJJ* 0,150006 01 





8LSD* 

0.166656 03 AVTEM* 0.198516 04 PH20* 0.468866101 

SSF* 

. 0,936496.03 

EMIS* 

0.429126-00 


QRAO* 0.417296-03* QRADS* Q.31536F-02 

J* 12 

XDtJJ* 0,160006 01 





BLED* 

0.166656 n3 AVTEM* 0.198276 04 PW20* 0.467796-01 

SSF* 

0.106966.02 

PM IS* 

tT. 4291 46*00 


— QRAn* 0.420256-03 QRADS* 0.357386-02 

J* 13 

XDtJJ* 0.170006 01 





BLED* 

0.166656 03 AVTEM* 0.198086 04 pw2o* 0.46*896101 

SSF* 

0.120356. 02 

PHIS* 

0.429156-00 


QP A n * 0.42102E-03 ORAOS« 0.39949F-02 

J* 14 

X0[J}| 0 # tPnonB 01 


FIGURE 11. RADIATION PROGRAM, SAMPLE PRINT OUT (Continued) 
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8LFD* 0.16665? 03 AVTgM* 0.19793E 04 pugn* 0.466l5g:01 SSF* U. 133746*02 PM I S* 0.42916F-00 

QBA>i« 0.42008E-03' QPAD3* 0.44149F-02 

J* 15 

xd { j i « o,i9nnne 01 

BLFD* 0*16665? 03 AVTE** 0.19781? 04 pP2n* 0.46556E-01 SSF* U. 14710?. 02 PM IS* 0.429l5?-00 

qqAo* 0.41776E-03 Q W A DS S 0.48327F-02 

J* 16 

XD t J ) * 0.20nnnE 01 

BLED* 0,16665? n3 AVTEM* 0.19772E 04 PM20* 0.46507?. 01 SSF* l< . 16037?. 02 FM I S* 0.4291 5F-00 

QPAu* n,41438E-03 QPAOS* 0.5247lp-02 

Ja 17 

xD ( j i ■ n. among 01 


SLED* 0.16665F 03 AVTEMa 0.19765E 04 pW2fl« 0.46470?. 01 

QP A U , 4 1 0 2 0 E - 0 3 QPAn$= 0.56573F-02 

Ja 18 

XOIJJ* 0.?2nnnE 01 


SSF * 0.17353F-02 PMIS- 0.42914F-OQ 


BLED* 0,16665c n3 AVTEMa 0*19759? 04 pM2n» 0.46441?^01 SSF* U.18654F-02 FMIS* 0.429l4g«00 

qpAp* 0.40534E-03 ORADS* 0.60626F-02 

Ja 19 

XDtJJ* 0 .230006 01 

*tFD» H , 1 6665 e 03 AVTEM« 0*19756? 04 PM2n» 0.46420E-01 SSF* . 19939F-02 FMIS* 0.42913F-00 

OR A *ia 0. 399986-03 QRAD$= 0.64626**02 

Ja 20 

X D ( J 1 ■ O.24P0OE 01 

BL E D* 0, 1 666 5? n 3 A-V^6** frHr97*3g-~ 04- -0-, 464 05E.01 SSF*--+i , 22397F-02 ?MI$a 0.4291 3E-0 0 

ORAna n , 76458E- 03 QP A OS* 0.72272P-02 

J* ?1 

XDtJJa 0.R6P00E 01 

BLFO* U.T6665C r)3 AVTEMa 0.19748? 04 PN2n* 0.46376?:01 SSF* U.24783F-02 FMIS* 0.4291 2F-00 

OPAMa 0*741496-03 QR A DS» 0.79667P*02 

Ja 22 

XD [ J 1 a 0.28nnnE 01 

BLED* 0.166656 03 AVTEMa 0.19744E 04 PMgn* 0.463506-01 SSF* U.27092E-02 FMIS* 0.42911F-00 

OBAO* 0.71673E-03 QRA0S= 0.86B54F-02 

Ja 23 

XDtJJ* 0 . 3 n n n o E 01 

9LFD* 0.16665C n3 AVTEMa 0.19743E 04 pM2n« 0.46347?:ui SSF« 0.2*3206-02 FMI$a 0.429UE-00 

OoAn. n. 691596-03 QP ADS* 0.93770P-02 

J* 24 

X D ! J ) ■ O*32n0O6 01 

8LED* 0*16665? 03 AVTEMa 0.19745E 04 pM2n* 0.46357E-01 SSF* i),30422F-02 FMIS* 0.429H6-00 

0 P A 0 a n , 34223E-0 3 QRAOs* 0.97192F-02 

Ja 25 

XD 1 J 1 a 0.33PD06 01 


BIRD* 0.16665? n3 AVTEMa 0.19776E 04 pM2n* 0.46648F-01 SSF* 0.320276-02 FMIS* 0.42938F-00 

OP A ita 0*501736-03 OR A DS * 0.10221F-01 

a BY AVT ♦ A VP METHOD 


TOTAL oADtANT HEAT F|_Ux IS 0 . 1 0220942E-M RTU/SFC-SQ FT 


FIGURE 11. RADIATION PROGRAM, SAMPLE PRINT OUT (Continued) 
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CWM* ( X ) 


CHAIN (4t 

biockaoe correction factors 


p 


a 


Mfcar not sr arsorrt t 0* method 
j« i 

xd t J ) • 0.56101E On 

0*» D>« U.29731E-n5 1ROAX* 0 . 2«731f-n5 

J« ? 

XDt JJ* O.6O10OE on 

0 R D •’ « G.51711E-04 ORDAX* 0.5*6*58-84 

J« X 

xorju 0,701006 on 

QRDab 0 , 15404E-03 ORIA*. 0.2f1872F.*3 

J* * 

*DtJI» O.SOnofiE 00 

qpD«» 0,270098-03 ORDAx* 0.*79728-rt3 

J- * 

XDtJJ* 0,900006 on 

0RD** 0.473326-03 0«1AX« n.9530*8..ft3 

J* * 

*D I J t » 0,100006 01 

0.347106-03 oRDAx* 0.13001E-02 

J* 7 

xdiju n.imooe oi 

QOD»« 0.165096-03 ftROARp 0.146526-62 

J* ft 

*D«J1» 0.121'OOE 01 

OoU“» 0.152636-03 f)RnAx* n,i*i79F-ft2 

J» 9 

«D ( J > • 0 .131006 01 

0®tM« 0.144976-03 nROAq. <5,176208-02 

Jp in 

*DtJ1« 0.14noo6 01 

QRD'-p 0.139316-03 nROAx* 0.190216-02 

J" 11 

XD I J J • 0,150006 01 

OBDap 0. 134306-03 oRDAXp n.2o36*F-o2 

J« 17 

XDUIp 0.160006 01 

QbDa» 0.131166-03 o«r)M* 0,2l6768-n2 

Jp 13 

XBtJI* 0.170006 01 

Q9B'*« 0.127536-03 oROAx. 0.279516-02 

J» 14 

XD1JJ« O.lRngpE 01 

0»Da* _Q. 126766-03 oRDAr* 0.242196-02 

J" 15 

0.190006 01 

ORUi* 0.125366-03 oROAp. 0.254726-02 

J* lo 

XD I J J • 0, 200006 01 

QoDa* 0.123446-03 OROAX. 0.267076-02 

J* 17 

XDIJ}« 0.2inonE 01 

OBD*b 0.121736-03 ORCA** 0.779246-62 

Jp 10 

XDtJjp n. 220006 ot 

0®D*P 0.119046-03 oRO A*» 0.291156-02 

J* 19 

XDIJIp 0,230006 01 

0®Dtp 0.117036-03 ORDAX* 0,302058-02 

Jp 20 

xQtJ)« 0. 240006 01 

ORD»« 0.222466-03 flROAp. 0.325146-02 

J* 21 

XDfjt« 0,261006 01 

ObD<p 0.215116-03 oRnA*. 0.346658-02 

• J* 22 

xDIJtp 0,261006 01 

OPDap 0.208666-03 #}RDAp« 0 . 3*753E-n2 

JP 23 

xDCJIp 0.301006 01 

OBDap 0. 201368-03 o«OAp» 0.387678-n2 

JB 24 

XOtJJp 0.320006 01 

OADa* 0.99589E-04 rjRDAp. 0 . 39763f?-o2 

J- 24 

XD(J)p 0.331006 01 

OBDa» 0.145216-03 ORO*P« 0.412158-02 

mE*t flu* «y absorption method 

TOTAL RADIANT -HEaT FLUX IS 0 .<121 49811-n? RTU/S8C-SO rT 


LEFT SfOB OF Wcg T 0> 
XOp 0.32X9867OF 01 
MP'DtE OF REGION 
XOp 0.27U695P6 U1 

right sins of *Exr.N 
XOp 0.32*986706 01 
QMO» ■ . 

OR1AP -. 

LEFT SIOF of HCGtO 

XOP 0. 310190376 05 
MI101& OF P6PlOfc 
XOP Q. 219*94098 01 
R 1 <shT Pinp Of ftf^i n 
XO» 0.3lol9Pj7p oi 
OR&Dp J, 

OROA" II. 

LEFT stop nf HC6I0 

XOp 0.307669836 (j 1 
M!noL6 OF PEQInjg 
XOp 0. lXo30413g 01 

RIGHT Pine or ngRl’N 

XOp 0.307*69838 o’ 

OR » DP 1». 

3R3A* v. 

LEFT SIDE OF RFC T 0 
XOp 0.297«479?F Cl 
RlnQLE OF RE6 I ON 
XO* 0. 89i69564fe 

XOp 0.11765444.8 P» 
*n» 0.11 *854448 M 

XOp 0.117654458 0! 
XO« 0.1178544*8 0 ) 
XOp 0.1176544-8 01 
Rir.HT P ] np OF »8Rl"N 
XOP 0.29744797 F IU 
ORADp ( 1.4350711*8-04 
OROAp (!.e»)3.*156lE-04 

left sids of rfgio- 

XOb 0,2878114078 03 
*1 IDLE OF REGION 
XOp O.Fft-OfcOSPf-O 
XOp 0,64?4333f!8 0 
XOp 0, 649433308 0 > 
XOp 0.649433308 0 
XOp 0.649433316 0 
XOp 0.649*33316 0 
RiSHT RlOp OF -i6ft) h) 
XOp 0.287604078 Ot 
ORADp 0,10821*7*6-03 
OROA* 0.15777*2*8-03 


J» 1 

0.13755223B-02 

TOP -0 , 1327131 1F-05 

Tl« -0. 137552? 3P -02 
ORADS« 0. 

ord as* o. 

Jp 2 

TIP 0 , 90989337P-02 

TO» -0.104045*4P-05 

n» -0.93989337P-02 
ORADSp 0. 

QRDAS> 0. 

Jp 3 

Tn« 0,274210*7P-01 

TOP -0.6«176»*5f-06 

TOP -0.27421087P-01 
QRAOSp 0. 

CrRDiS* 0. 

Jp 4 

TOP 0 , 46014 039F* 01 

YO« -0.4364(54788-06 
TOP -0,576791*4P-06 
TOP -0.576791*48-06 
TOP -0.57679114F-06 
TOP -0.57679154P-06 
TOP -0.57679154P-06 

TOP •0.46014039P-01 

0RAOS- 0.835071138-04 
QRDA$b 0.9033l56lg-04 
Jp 5 

TOP 0 . 620 103*2P-01 

TO* *1.237887*38-06 
TO* -0.314413798-06 
TOP -0,314413298-06 
TOP -0.314413298-06 
TOP -0.3144l3?gp-06 
TOP -0.314413298-06 

TOP -0 . 620103*28-01 

ORADSp 0.191725888-03 
CRDaSp 0.24810*828-03 


FIGURE 11. 
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LEFT STDP OF R"G10 


J* 10 


J* 14 


XO* 0 .27*847040 01 
MIDDLE OF »EGIon 
XO* Q.154«251«F-0 
XO* 0.20x5>73446.0' 
XO* 0 .203573440-0 ' 
XO* 0 . 20*573440-0 ' 
XO* .0, 20*573440-0 
XO" 0,<? r| li > 7344 & _(r 

RIGHT SlOC OF RERl'N 
XU* 0.278*470«F 01 
f)R* D « ;i.l2379'i246-Q3 

QRn*» n. 1 1569977E-U3 

LEFT S U10 OF R0GIO 
XO* 0.27P50731E 01 
MIDDLE OF REGION 
XO* -0 .153176340-0'’ 
XO* -0.201 129R4F-0'' 
XO* -0.2nil2OH4F-0. 
X0* -0.2111129840.0 
XO* -0.201 129M4F-0" 
XO* -0.201 12944F-0 ■ 
RIGHT R I np OF RED I 'N 
XO* 0 , 27(180 731 F 01 
QR a 0* is , 1 22*0 92 it - 03 
ORDA* ’>.*5O29124fc-04 

LEFT SIDE OF HFGIG' 
X0« 0 . 26?7250l)F 0) 

MIDDLE OF REGION 
*°« -0 # 3<»42763lE.0i' 

XO* -0.521 11 OHOF 0" 
XO* - 0 . 521 1 1 OflrtE 0 ; i 
XG* -0 . 521 HORDE 0 1 
XO* -0.521 linanp 0 ■ 
XO* -0.521110800 0 ; , 
RIGHT Sine OF RFGl’N 
XO* 0 . 2637250 OF 01 
QP*0» 12944^7*6-03 

QRDA* I! .MU8762440-O4 

LEFT STOP OF RPGTO'' 
XO* U.25D45007fc 01 
Mlnoi.E OF REGMn 
XO* -0.6031446RF 0 1 ' 
XO* -0.79716991F Or 
XO* -0 . 79716991 F O' 
XO* -0. 79716991 F Q'> 
XO* -0 . 79716991 t 0 " 
XO* -0 . 7971 6991 f O' 
RIGHT SIDF OF •> E G I " N 
XO* 0.25545007F 01 
QR*D* n ,134096156-03 
QftDA* n.4«3?2777t-04 


TO* 0.759303480-01 

TO* -0. 753814*50-07 
YD* -0 .996308" OP-0 7 
TO* -t) . 996308*00-0 7 
YD* -rj.99f*3Q0POP-f)7 
To* -0 ■ 996308"0P-'J 7 

Y0 * -0 t 99630 rh0 f .U7 

YD* -0 . 759303480-01 

QRADS* 0.315516110-01 
QRDA5* 0 .3*380459e-03 
J* 7 

TO* O.R0363512P-Q1 

TO* 0, 744766230-07 
YD* 0.98434949F-07 
YD* 0 . 984349M9F-07 
TO* n , 98434959F-07 
To* 0,964349590-07 
YD* 0 , 9843 4 9"9F -G 7 

TO* -n. 8838351 2F-01 

ORADS* 0,438125340-01 
QRDAS* 0.41883 3 91P-01 
J« 8 

YD* 0.996302070-01 

Yn * 0,19296277^-06 
YD. 0.255036«4P-l'6 
YD* n.25%036«4P-06 
YD* 0.255036 8 4F-06 
YD* 0 . 255036" 4F-06 
YD* 0.255036*40-06 

YD* -0 , 9963020 7P- (11 

ORADS* 0.58 777211 F- 0 1 
QRDAS* 0.4*9710160-01 
J* 9 

YD* (1.109994*50-00 

YD* 0.29518505F-06 
YD* D.39014294P-06 
YD* 0. 39014294F-06 
YD* 0.390142Q4F-06 
YD* 0.39014294F-06 
YD* D.390142O4F-06 

YD* -0 » 109994*50-00 . 
QRADS* 0 . 7niA6fl26F-fl3 
QRDAS* 0 .5180329*0-1)3 


LEFT StUG OF RRGTO'- 
XO* 0.24883580E 01 
MIDDLE OF REGION 
XO* -0, '79?7869 3 E 0 
XO* -0.10478180E 01 
XO* -0.1047818DF 01 
XO* -0.1 04781806 01. 
XO* -0,104781806 01 
XO* -0.104791806 01 
RIGHT SlDF OF REG I f.'AI 

X°* * E 

0 2406358® 01 

QR A D* 0,1 371861 4fc-03 
QR D A* 0,4643635 3 fc-04 

LEFT SIO0 OF REG TO-'- 
XO* 0. 24??4354E 01 
MIDDLE OF REGIDN 
XO* -0 , 96*31 82?F 0° 
XO* -0,127^81596 01 
XO* • 0 . 127^81596 01 
XO* -0.127981596 01 
XO* -0.127981596 01 
XO* -0.127981596 OJ 
RIGHT SIDF OF aEGl'N 
XO* 0.24224355E 01 
QHAD* 11.1 39097876-03 
(3RD** '1.447660316-04 

LEFT STO0 OF RPGIU- 
XO* 0.21623776E 01 
MIDDLE OF REGION 
XO* -0.111934576 01 
XO* -0 . 1*7942726 01 
XO* -0.14794272E 01 
XO* -0.147942726 01 
XO* -0.14 7 94272F 01 
XO* -0.14 7 9427?F 01 
RIGHT RIDF OF RES I °N 
XO* 0.23623776E 01 
QP AD* U. 1*0084546-03 
QRDA* 11.43719(1806-04 

left SIDE OF RPGIO^ 
XO* 0 , 23n 58719F 01 
MIDDLE OF REGION 
XO- -0,125297298 01 
XO* -0.165604096 01 
XI • -0 , 16460409F 01 
XO* -0.16*604090 01 
XO* -0.16560409F 01 
XO* -0.165604090 01 
RIGHT 8lDP OF ->EGl N 
XO* 0,210487196 01 
QR»D* 0. 1403*1496-03 
OHDA* J.42508919E-04 


YD* 0.11956616F-00 

YD* Q.387997" 60-^6 
YD* 0.512812646-06 
YD* 0.512812"46-06 
Yd* 0.512812640-06 
YD* 0.512812640-06 
YD* 0.512812640-06 

Y n- ' ■ 6- 

o 119566-76 00 

QRADS* 0 .819054400-01 
QRDAS* 0,564469260-01 
J* 11 

YD* 0.1 28493950- 0 0 

YD* Q.47390463F-06 
YD* 0 .626354700-06 
YD* 0,6263547(30-06 
YD* 0.62635470P-06 
YD* 0,626354700-06 
YD* 0.626354700-06 

TO* -0*12849 3 95F-00 

QRADS* 0.970152270-01 
QRDAS* 0 . 6 n 9?3532P- 0 1 
J* 12 

YD* 0.13686175C-00 

YD* 0.5*7818970-06 
YD* 0.724046400-06 
YD* 0.72404640P-06 
YD* 0.724046400-06 
YD* 0.724046500-06 
YD* 0.7240*6506-06 

YD* -0.136861750-00 

QRADS* 0.111021680-0? 
QRDAS* 0.6*295*340-01 
J* 13 

YD* 0.144747060-00 

YD* n .611217*00-06 
YD* 0.810483 1 ' 10-06 
YD* 0.810483010-06 
YD* 0.81 04R3O10-06 
YD* 0,810483010-06 
YD* (1.810483010-06 

YD* -0.144747060-00 

QRADS* 0. 135657810-03 
QRDAS* 0.695463310-01 


LEFT Strip OF RCGIO . 
XO* 0.225263910 01 
MIDDLE OF REGION 
XO* -0.1 3 7 091 2lF 01 
XO* -0.1 *1 191990 0 1 
XO* -0.101191990 01 
XO* -0.101191990 01 
XU* -0,101191990 01 
XO* -0,181191990 01 
RIGHT RIDF OF REG I'M 
XO* 0 , 22526391 0 01 
QFUG* 0,1 4002769E-03 
QRDA* 0,422596846-04 

LEFT S f Q0 OF R0GIO- 
XO* 0,22 o 2423D6 01 
MIDDLE OF REGION 
XO* -0.1*7959200 01 
XU* -0,19*556226 01 
XO* -0,19*556226 01 
XO* -0,19*556236 Ol 
XO* -0, 19*55622E 01 
XO* -0.19*556226 01 
RIGHT «IDF OF WEGI^N 
XO* 0.22n2423ng 01 
OR AD* 1.i 39254656-01 
OHDA* U. 4178535(16-04 

LEFT SIDE OF R0GtQi 
XO* 0.21*5006.16 01 
MIDOLE of region 
XO* -0.159917806 01 
XO* -0 , 21 j 36165E 01 
XO* -0.2H 361656 01 
XO* -0 , 21l 361656 01 
XO* -0.2113616*6 01 
XO* -0 .211 36165F 01 
RIGHT *ID0 OF REG l iN 
XU* 0.21*500636 Ol 
QHAD* ", 1381281*6-03 
DRDA* ■' . 41147?66t-04 

LEFT ST06 OF RFGtOJ 
XO* 0.211019136 01 
MIDOLE OF REGIdn 
XU* -0.1659326*6 01 
XO* -0.21931146F 01 
XO* -0.21931146F 01 
XO* -0.219311466 01 

XO* -0.210311460 oi 

XO* -0 . 2 lo31 14*F 01 
RIGHT Slop OF «FGI'N 
XO* 0.211019126 01 
QHAD* ".1 36731976-03 
QROA* i .405771*16-04 


YD* 0.152213130-00 

YD* 0.670938100-06 
YD* 0.88677174P-06 
YD* 0.886771740-06 
Yd* 0.886771740-06 
YD* 0.886771740-06 
Yo* 0.886771740-06 

Yd* -0.152213130-00 

QRADS* 0 .119860610-0? 
QRDAS* 0,717723000-03 
J* 15 

YD* 0.159303930-00 

YD* 0,724127*00-06 
YD* 0.957071730-06 
YD* 0.957071730-06 
YD* 0.957071730-06 
YD* 0.957071730-06 
YD* 0.957071730-06 

YD* -0,159303930-00 

QRADS* 0.153786066-03 
QRDAS* 0 .7795O8250-D3 
J* 16 

YD* 0.166078360-00 

YO* 0,782653730-06 
YD* 0.103442510-05 
YD* 0.103442510*05 
YD* 0.103442*10-05 
Yn* 0.10344251F-O5 
YD* 0.103442510-05 

YD* -0.166078360-00 

QRADS* 0 . 167598880-07 
ORDAS* 0 .82065551F-03 
J* 17 

YO* 0.172583020-00 

YO* 0.81209117P-06 
YD* 0.107333230-05 
YD* 0.107333230-05 
YO* n. 107333230-05 
VO* 0.107333330-05 

Yn* 0.107333330-05 

YO* -0. 1725830 2f -00 

QRADS* 0.18137204F-03 
QRDAS* O.8MP32670-O1 
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J* 18 

left sioib of BFsto • 

XO* 0.20*7782*6 01 YO* 0 . 1 788331 4F- 00 

hipdie or REGinn 

XO* -0.173737926 01 YO* 0 .650290788-06 
XO* -0.2296275*6 01 YO* 0 . 11 2382036-05 
XO* -0.279627566 01 Yn* 0.112362P36-05 
XO* -0.2296275*6 01 Y<)* 0.112382536-1)5 

XO* -0.2296275*6 01 Vn* 0.112352136-05 
XO* -0.229*275*6 01 Yn* 0.U236213F-05 
HlfiWT S|np OF *»et3 1 >N 

XO* 0.28*7782*6 0J YO* *0 . 1 788331 46-00 
<3M0* ». 135112296-03 0»Af>8* 0 . l«47833lf-0? 

DBDA* 1.396790946-04 QRDAS* 0 . 90 091l7*p-03 

J* 19 

UFT STOP OF »6Gt0‘ 

XO* 0. 20?76?056 01 TO* II . 184868496-00 
Hjnote or PtGIOfc 

XO* -0 . 18fl9382**6 01 YO* 0.885530056-06 
XO* -0.23ol441«6 01 YO* 0.117039566-05 
XO* -0.23O144196 01 VO* 0.117Q3956P-05 
XO* -0.230144196 01 YO* 0, 117019566-05 
XO* -0 , 23ol441«F 01 YO* 0.11703956F-O5 
XO* -0.239144196 01 YO* f> . 1 170395*6-1)5 
BIGHT 5 1 OP OF "MUN 

XO* 0.207762056 01 YO* -0. 184868496-00 
(}8*D* <.33332*7*E-03 QBAOs* 0 . 2*61 1594C-03 

0«PA* ii. 39 0089576-04 QHD*S* 0 . 93992l)7?p-l»3 
J* 20 

LEFT S ! 06 OF RCfitO 

XO* 0.19*955786 PI YO* 0 , 190720 (S9F-00 
HjnotE Or RECJOt. 

XO* -0 , 1 54*18906 01 YO* 0,914310416-06 
Xfl* -0.24*916536 01 TO* 0.120843436-05 
XO* -0.24*916356 01 YO* 0.120841436-05 
XO* -0.246916556 01 YO* 0 . 1 20#43«36-05 
XO* -0.Z4691655F Cl YO* 0.120843436-05 
XO* -0 .246916556 01 YO* 0.120843436-05 
RIGHT 3 1 06 OF »6«l 'N 

XO* 0,19*9557*6 01 YO* -0.190720*96-00 
0940* 1.954*5*936-03 08*03- 0.233*018*6-09 

ORTH* (1.7429189*6-04 0*0*8* 0 .1*14212*6-03 

•J* 21 

L6FT Sf 06 OF «6GJ0- 

XO* 0 . 1 *1 922026 01 YO* 0 •'202000*16-00 

hi note of »eeio*« 

XO* -0. 197*60266 01 YO* 0.967369116-06 
XO* -0.261245466 Ol YO* 0.127856146-05 
XO* -0,26i2454*F 01 YO* 0.127856146-05 
X0« -0.2*12454*6 01 YO* 0,127856146-05 
XO* -0.2612454*6 «1 YO* 0 . 1278561 4F-05 
XO* -0.2612454*6 !U YO* 0.127856146-05 
BIGHT * I rip or «6ftl‘ N 

XO* 0 * 19i 92202F 01 Yi* -0.202000*16-00 
<JR*0« ». 347169846-03 QBAOS* 0 .25fl«*l«,e-fi9 

0604- ! ,717047796-04 Q*0«S« 0.1118591746-03 


U6FT STOP PF B6GT0- 

X0» 0.15857304F 01 Yn* 0.212729*76-00 

HJPD 16 or H6GIOH 

XO* -0.20*161586 01 YO* 0.100897506-05 
XO* -0,273*814*6 01 YO* 0.133355156-05 
XO* -0.2735814*6 01 YO* 0.133355156-05 
XO* -0.273481**6 01 Yn* 0.133355156-05 
XP* -0,2734614*6 01 Yn* 0.133355156-05 
XO* -0.2734814*6 01 Yn* 0.133355156-05 
BIGHT 51*16 or 96Gl«»* 

XO* 0.15*573046 01 TO* -0.212729376-00 

GRID* <.338910656-03 ORAOS* 0 . 2»27Q9256-P7 
<jRt,** 'I. *96210276*04 Q«D*S* 0 . 11555384C-0? 

J* 23 

LEFT ST06 OF H*GtO 

XO* 0.179821656 01 Yn* 0 .2229711 76-00 

HlPOtfe OF HgGinN 

XO* -0 , 213 1679 *6 01 YO* 0.104277*66-05 
XO* - 0 . 2 * 161 ) 95*6 Ci n« n. 137522556-05 
XO* -0.2*1 * <195*6 01 YO* 0.137822*56-05 
XO* -0.2*1*f*95 B € PI YO* 0.137822*56-05 
XO* -0.2*1*095-6 01 Yn* 0 . 137622596-15 
XO* -0.2*1*095-6 Oi TO* P . 1 37B?2*56- i, 5 
BIGHT Bine of i*60} n 

XO* 0.179*217*6 01 Yn* -0.222971176-00 

gp»0* -.?305294Pfe-03 QBAOS* 0 . 30576 219P-09 

IRn** :.i .*71214406-04 QRDAS* 0.12226899P-03 

J* 24 

UfeFt S 1 06 OF B6G»0 

XO* 0.174*98956 Ci Yo* 0.232433576-00 

Hi note 

XO* -0.217*52556 <M Yn* n .106619226-05 
XO* -0.2*7933396 01 to* 0.140917496-05 
XO* -0.2*7933396 01 Yn* 0.140917496-05 
XO* -0.2*7933596 01 Vo* 0.140917-96-05 
XO* -0.2*7933336 03 Y«* 0.140917496-05 

XO* -0.2*7933396 01 YG* 0.140917496-05 
BIGHT 5 1 G6 OF 06G1--H 

Xt)» 0.174*9*9*6 01 Y n* -P. 232833*76-00 
OHIO* U. 11407*376-03 0B408* 0.31666973P-83 

GRPA* P. 331964306-04 QB0»S* 0.125565636-03 
J* 25 

terT STOP PF BFG10 

XO* 0.1731672*6 01 Yn* 0 .23765675P-00 

H I HQL6 or B 6G l 

XO* -0.220407256 01. YG* 0 . 1 0 78695 1 C-05 
XU* -0 .29i 3(;9916 01 YO* 0.1*2559996-05 
XO* -0.291309916 01 Yn* *.l4?5*909e-P5 
XO* -0 . 291 309916 01 TP* 0.142569996*05 
XO* -0 .291 30991 F 01 YO* fl. 1425*9996-15 
XO* -0.291309916 01 Yn* 0.142569996-05 
RIGHT «IOr OF WfGlH 

XU* 0.1721672*6 01 YO* *n . 23765675B-00 

fjfi »D« 1 . 1 672*371 6-03 OB*DS* 0.333394006-87 


asn** , i,4H4 02*5°fc-04 QRDAS* 0.1.30425656-37 

H1TH hi OOx»GF TPBB. CTInH.TME R Al) I AT | v6 Fi.iiv *Y TH6 AVT. »M) AYR. MFTMOO IS 0.333394 098-02 9TU/S6C-SO FT 
BY THfc a 0 onfiP T T O' 1 €TnnD» TT IS * . 1 30 42M86-G2 
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CM* INJ {?] 


M 0Y AVT + A VP HtTH'in 

rK*-o.i973^p m 
PhJ* n.4°Vspc*i>(j 

j= 1 

xnt ji * n .smme on 


PLCU* 

i .1 fl6(!6F.rU 

AtfTEM* 0 .33161 E 04 PH?ri = 

0.12«65f-(i0 

SSF* 

• 254 03F-05 

FMISa 

0.29740F-02 


fi»A 

* .43 4 986-06 G*Ai)Sa 0, 

j* 2 

xnrjls 0. A O' 1 one 

,43490p-nft 

on 





Rt>D* 


AvTfcM* 0.33569F 04 PM?ti 

0.13117E-U0 

SSF* 

".47398F-04 

FMIS* 

0.12041F-P1 


(;r> A 

* 1 . 29743 E -04 quads* n, 

J* 3 

X U ( J 1 a " . 7 1 1 M F 

,30176^-04 

00 





RIFD* 

,n229di- nii 

AWTfcM* 0 . J 6 4 7 *i F Q4 PM?n* 

0 • 1 4P65F - 0 

SSF * 

1.16325F-03 

FMISa 

0.229O3F-01 



* !. 239096-03 QRAiK* 11, 

J* 4 

X U 1 J 1 ■ n . * o *■ n o F 

.26227F-03 

00 





RLFO« 

.1 ?8R2 C "1 

AVTEm. 0.33750*? 04 pw?ns 

0.131 75F-lifI 

SSF* 

.36853F-03 

FMISa 

0.363M7F-nl 


fJOA - 

* '.43A60E-03 QwATSa 0, 

j* b 

xorjiB o.->nM*e 

. 70UMC-03 
on 





0L p D * 

i,l 749.tt "1 

AVTEM* f).304ft?F O'* P*- ? 0 s 

0.H144F-C0 

SSF* 

» .61387F-03 

FMISa 

0.4«9n3F-0l 


r;oA 

>■ '.*317236-03 chads* n t 

Js to 

VD ( J 1 a p , 1 ni.niF 

,12l4ir-n? 

(11 





ROD* 

!*.?4 n 4*c nj 

AVTE^a 0. 2*161 F (14 

0 * 960926 - (i 1 

SSF* 

: . 91699C-03 

PMIS* 

0 .633R5F-01 


CP* 

= .57629E-03 QkAdS* 0, 

J- 7 

X i.M J ) a 0.11 ' n o t 

, 1 ?tt V4C -0 ? 
01 





R l F 0 * 

n,9*t*D4F Cl 

4 V T6 M a 0 . ?6*1*F O 4 PHgoa 

0.«flo34F.ul 

SSF « 

■ * . 15588P-02 

FMIS* 

0. 7P5*!8F-nl 


CPA 

= i.6ln61E-03 (J^AOS^ * , 

Ja 8 

xntjja o . i 2 r n n f 

,?40noc-02 

01 





Rlf=0* 

•:.336«of m 

4VrtM« O.25203E 04 pi-?n* 

n.79473F-0l 

SSF * 

.16320&-02 

FMISa 

0 .an7l 0F-01 


Op A 

* '.S6595E-03 Op A y <5 r n, 

J* 9 

XDIJI* P .M , '0O f ; 

,?96ftOP“fl2 

01 





0LCD- 

■ i . 4 n699 c 01 

AVTgMa 0.24Q4OF 04 PM?n* 

0 . 71O02F-01 

SSF* 

i.2^306P-02 

FMISa 

0 . 9f 990F -0 1 


(j p A ■ 

* J.57715E-03 O*A0s= 0, 

J* 10 

XDIJ1* 0,1. An ME 

,3563lP-0? 

01 





Rl>0* 

■I.4S97XP «tl 

AVT6*« 0.?305f'F 04 pM?n* 

0.h5A23F-01 

SSF* 

’ .244966-02 

FMIS* 

0.973O2F-01 


go A 

= i .548496-03 QUADS* 0, 

Ja 11 

xDIJIi n . 1 5 f- o n E 

.4H16F-02 

01 





Rl c D ■ 

; 1 . 5 3 4 7 7 f r 1 

AvTEMa 0.22233F 04 pm? n * 

0.60693f:(;1 

SSF* 

.24850F-02 

FMIS* 

0.1O5O0F-00 


go A 

* ''.536336-03 Q»AD$s 0, 

J» 12 

XntJ)« 0 . 16 MnF 

, 46479F ” 0 2 
01 





RLFO- 


A J T 6 M * n . ?173^F 04 PMpna 

n .57094F-U1 

SSF* 

.3U31F-02 

FMIS* 

0.111*9F-00 


go* * r '. 532746-03 0 PA)Ss f .51 8"7c-02 
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J* 10 



XDUla 0 . i7nonE 01 





- 

0LED» O.652O0F 01 AVTEMa 0.212226 Q4 PH20* O.54998E-01 

QPAD* 0.519346-03 QRADS* 0.570006-02 

J» 14 

XD u i * o.tanonE 01 

SSF* 

0.379106-02 

6M!S* 

0.117056-00 

* 

BL60» 0.736086 01 AVTEM* 0*206996 04 PH2o* 0. 519006-01 

Q 9 A D ■ 0.506656-03 0RAHS* 0.620666-02 

J* 15 

XD u ) » o.i9non6 01 

SSF* 

0.425606-02 

PMIS* 

0.124666-00 


~9Lf=D* 0.001546 0l AVTEM* 0.205456 04 PH20* 0.510586-01 

QRAD* 0.52011E-03 QRADS* 0. 672686-02 

J. 16 

XD I J 1 ■ 0. 200006 01 

SSF* 

0.472596-02 

6MIS* 

0.130476-00 


9L6D» 0 . 092736 0l AVTEMa 0.202206 04 PWjrj* Q. 491056^01 

Q P A 0 • 0.51740E-03 QRADS* 0.7244JF-02 

J* 17 

— . XD ( J la 0.210006 01 

SSF* 

0.519916*02 

6MIS* 

0.137S9E-00 


BLBD* 0.987966 0l AVTEM* 0.20077E 04 PH2n* 0.482846-01 

QPAD* 0 .531006*03 QRADS* 0.77752F-02 

J* IB 

XDU1* 0.2200DE 01 

SSF* 

0.567366-02 

6M!Sa 

0.144066-00 


SLED* 0.108736 02 AVTEMa 0.19965E 04 PH2n« 0.476316-01 

QPADa 0. 544866-03 QRADS* 0.832006-02 

j. 19 

XDtJU 0,230006 01 

SSF« 

0.614816-02 

FMISa 

0.151786-00 

w 

8L6D* 0.114166 02 A VT£M» 0.1989P6 04 PM2n* 0.47?426-Gl 

QR A n» 0.549356*03 QRADS* Q. 886946-02 

J. 20 

XDUla 0.240006 01 

SSF* 

0.662176*02 

6MIS* 

0.155436-00 


BLED* 0.124846 02 AVTEMa 0.1984RE 04 pM2o» 0.469506-01 

QRADa (1.110456*02 QRADS* 0.997396-02 

J. 21 

XD t J ] ■ 0.76n00E 01 

SSF* 

n. 754106*02 

PHIS* 

0.162006-00 


BLED* 0.139636 02 AVTEM* 0.19775E 04 PM2n* 0.465256-01 

ORAD* 0 .113566-02 QRADS* 0.111106-01 

- J- 22 

XDUla 0,200006 01 

SSF* 

0.64465E.O2 

EMISa 

0.171076-00 


BLED* 0.152716 n2 AVTEMa 0.19722E 04 PM2n« 0.402186-01 

QRAD* 0.115266*02 QRADS* 0,122626-01 

J* 23 

XD ( J ! * 0,300006 01 

SSF* 

0.933936*02 

PM I S» 

0.179026-00 


BLepa 0.166656 n2 AVTEMa 0*197106 04 PW2n* O.40195E-O1 

QPAD* 0.117596*02 QRADS* 0.134386-01 

J* 24 

XDIJl* 0.320006 01 

SSF* 

0.102106*01 

6MIS* 

0.187076-00 

* 

"BLED 1 * (T.f 8t5<rr 02 AVTEM* 0.19737^ (74 PW20* 0.463076-01 

QPAD* 0.613646-03 QRADS* 0.140526-01 

J* 25 

X0U1* 0.330006 01 

SSF* 

0.106446*01 

IMIS* 

0.195036-00 

0 

8L6D*“lT.T8n<J6 02 AVTEMa 0.200136 04 P» 2n« 0.489226-01 

QPAO* 0.98517E-03 QRADS* 0.15037F-01 

Q BY AVT ♦ AVP m6THdd 

-WAL HEAT FUK 0 . 15036862E-01 8TU'SPC-SQ *T 

SSF - 

0,112806*01 

6MIS* 

0.202076-00 
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CHAIN U1 


meat flux er absorption method 

J* 1 

XD I J 1 ■ 0.56^006 on 

QBU« 0 , 746*95-07 oROABi 0 . ?4ft89F-fl 7 
J* 7 

xd i j ] * n.Aonoofe on 

QPD.ft* 0.50860E-05 oRDARb 0.51607E-05 
Jb 3 

XD I J 1 * 0.70000E 00 

0 p L> 4 ■ 0 , 36463E W 04 oRDABb n.4i623F.Q4 
J«i 4 

xdiju n.ROPone oo 

ORD/« 0.64390E-04 ORnARB 0 . 1 0601F-0 3 

Jb •> 

xd i j i • o . 90 n n oE on 

QBDab 0.80574F-04 oROABb 0.1R659F-P3 

J* 6 

XD I J ) • O.lOnooE 01 

QBDab O.93671F-04 (JROARb 0.2*02*F-n3 

Jb 7 

xD I jl ■ O.ltnoBE 01 

QBDab 0.95493E-04 ' OROABb ,0 . 375?5F«.f>3 

Jb 8 

XD I J I ■ 0 , 12 P 0 OE 01 

QBD * ■ 0 . 94769E-04 oRPARb 0.47052F-o3 

J* 9 

XDU1» 0,13nonE 01 

QBDab 0.95205E-04 oRPARb 0.5ft573F.p3 

Jb m 

XD ( J ) ■ 0,14000E 01 

0»D:.« 0.92091F-04 oRDARb 0 . 6R782F-fl3 

J« 11 

XD(J)« 0.19000E 01 

QR0a« Q.91110E-04 n«0A9* 0.74693F-p3 

J* 12 

xd( J)* o . 1 6 n o n E oi 

QBDab 0.886IS8E-04 oROABb 0.83759F-P3 

, J« 13 

*D * J 1 * 0.17000E 01 

ORD 0.86632F-04 f)ROARi 0.9;>623F-n3 

J* 14 

xd ( J i ■ o.ienonE oi 

ORDAt 0.90721E-04 oHP AB b 0.ini69F«n2 

JB 15 

XD(J)b 0,19f>0fjE 01 

ORDftt 0 . 91227E- 04 oRPASb O.11O02F-P2 

Jb 1 A 

X D ( J I > 0.20000E 01 

OBD*« 0.94834E-04 oRDARb 0.12030F-P2 

J» 17 

XDIJIb 0 , 21 0 0 f)E 01 

GBDab 0.10179E-03 • oRDARb 0.13048F-Q2 

Jb 18 

XD I Jl * 0.22000E 01 

QBDab 0.11091F-03 ORDARb 0.14157F-n2 

J- 19 

XDIJJ" 0.23000E 01 

OBD ab 0 , 11596E- 03 QRDARb 0.153176-02 

Jb 20 

XDIJlB 0 , 24 n ooe 01 

QBDab 0.25189E-03 oRPARb 0.l7836F-f»2 

JB 21 

xorjji o ,2fenonE oi 

QBDab 0.28888E-O3 ORPARb 0.2fl724F-p2 

JB 27 

XDIJ1> 0 • 28(1 0 oE 01 

ORDAb 0.32669E-03 ORDARb (1.2399lF-o2 

JB 23 

XDUIb 0 .300 DOE 01 

QBDab 0.37396E-03 ORDARb 0.27731F-02 

J* 24 

XDIJIb O.32O00E 01 

QBDab 0.22030F-03 ORPARb 0.2»934g-n2 

J- 28 

XDIJIb 0.33000E 01 

QBDab Q.55214E-03 oROArb P . 33455F-H2 

meat flux ay absorption method 
TOTAL BADtANT MEAT FLUX IS 0 . 33495?45E-0? RTU/SFC-SQ FT 
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CHAIN ( * 1 

fcfiOC*»GE CORRECTION FACTORS 




r 


* 




LEFT side of rpgto 
XO* 0./5rh651«F 01 

middle of region 

XO* 0.75R51804E 01 
SIGHT Slop OF wpGl'iN 
XO* 0./5oll755E 01 
OHIO* il. 

0*0*. .1, 

LEFT S1DF OF HCGIO 
XO* 0.74aV 1Q21E 01 
H I 001 E OF KESION 
XO* 0. 74<t000*F 01 
SIGHT sine OF "»FG I N 
XO* 0.7485776OE oi 
QR A0* <1, 

r>sn *• i.i, 

left STUB OF SBGIO 
X0" 0.M8A2355F 01 

HIODLE OF HER I ON 
XO. 0.7ti5596?E 01 
SIGHT Sine OF «EGl-’N 
XO* 0.72 .g6322*E 01 
QR«D> R • 

QRDA* 0. 

LEFT 5 f OF OF S B G 1 0 ' 
XO. 0. 69*318416 01 
MIDDLE OF HER I ON 
XO* 0.68 fl 2874SE 01 
SIGHT Sine OF ^EGl ’ N 
xo. 0./ftft.913'’E 01 
OS t0> ii. 

QRDA« 0. 

LEFT SlOF nF RRGTO 
XO* 0.66.00777E 01 
H I DOLE OF HER I ON 
XO. 0,65i58.7.E 01 
RIGHT sine OF HERl'N 
XO" 0.67o2564«E Cl 
ORtQ. 0. 

0*0*. il . 

LEFT SIDE OF SFGtOv 
XO* 0.64544933E 01 
SI COLE OF REGION 
XO. 0. 62511*5. E 01 
SIGHT SIOF OF WERl'N 
XO. 0 . 65ol 867 nE 01 

OS tD* ‘J. 

osn*. n. 


J* 1 

Tn. 0,202822324 01 

TO. O.20222274F HI 

TO. 0.20187571F 01 
QRADS* 0. 

osoas. o. 

J. 2 

TO. ft. 207480.4* 01 

TO. ft . 1988121 SB 01 

TO* ft • 1 9. 2 2639 B 01 
OR AOS" 0. 

QSDAS. 0. 

J. 3 

TO. 0.201677A5F 01 

TO. 0.189703*2F 01 

TO. 0.182892984 01 
ORAOS* 0. 

QRDAS. 0. 

J> 4 

TO* 0.2Q147374P ftl 

TO* 0.181366394 ill 

TO. 0 . 1 700827SP 01 
ORAOS* 0. 

QRDAS* 0. 

J* 5 

TO. 0.20ft 83375F 01 

TO* ft. 173714. ie ftl 

TO* 0 . 1 58679. OF ftl 
QRADS* 0. 

QSDAS* 0. 

J. 6 

TO* ft . 1998891 4C ftl 

TO* 0,16665759* ftl 

TO. 0.14843652* 01 
QR AOS* 0. 

QRD*S» 0 • 


LETT STOP OF SPG 10 

XO. 0.62443291E PI 
MIDDLE OF REGION 
XO. 0 ,6ftn5335?E 01 
BIGHT Slop OF HEfil .* 
XO> 0.64n38663E 01 
QR/iQ. ft. 

QRDA* .i, 

LEFT ST’JF OF RECTO 
XO. Q , 60*7978. E ftl 


J* 7 

TO. H.19881179B 01 

TO. 0. 160103. 9F 01 

TO* 0.139095S0P ftl 
GRADS* 0. 

U»DA5« 0. 

J* 8 

To* 0.1976601 4B 01 


HIDDLE OF nEG t dr 
<0* 0,5776«M37E 01 

SIGHT Sine OF HFRl-'N 
XO* 0.62974721* 01 
QStQ. 

OR DA* u. 

<*§fT Sfpe n r S|»G| 0 
XO* 0.5864G474E 01 
HI DOLE OF H§GI"N 
XO* 0.55.13224E 01 
SIGHT sine OF >»ER I f1 N 
XO* 0 , 60.180 7*E 01 
QR t Q. ft. 

USD A* ft. 

left SIDE OF RFSTO- 
XO. 0.?6ol3869E 01 
MIDDLE OF R6Gl"N 
XO. 0.53S97952E 01 

right sine or *pRi f 'N 

<0. 0.59&59S41E 01 

QRaO. a. 

QRDM 3. 


T". n,lS399ft«9F ftl 

To. it . I3ft921.se 01 
QUADS* 0. 

OSD AS* 0. 

J. 9 

TO. 0.196499P7F 01 

TO* 0.148266204 01 

TO. 0.12218724* 01 
QRADS* 0. 

QROA$* 0. 

J* 10 

TO. 0.19S34531F 01 

TO* 0.14269344F 01 

TO* 0 . 1 152201 IF 01 
QRADS* 0. 

QRDAS* 0% 


J« 11 

T 8 OH 0 

LEF IDF F *GI 

XO. 0.55*ft9?3?E ftl TO. 0.194J2529F 111 
MIDDLE OF REGION 

XO. 0.M7O103SE 01 To« 0.137#36?1P 01 
RIGHT Sine OF ft E 11 1 • N 

XO* 0 . 67*91 71. F 01 TO. ft.lO«J43.3F ftl 


ON *0* 0, 

QRADS* 0. 

GRDA* n. 

QRDAS* 0. 


J* 12 



XO. 0 . 5375931 **E 01 

Ti. ft . 1 9313853** 


H1O0I E OF HER I in 

XO* 0.49O11991E 01 TQ. ft. 133Q6658B Ml 
SIGHT Sine OF «ERl‘N 

<0. 0.56206989E 01 T<J. 0.KI1»1J*1P 01 

QR tD* ft. 0**ns. o. 

QRn A* l], QRDAS* 0. 


LEFT Sfpp Op Rkgi^ 
XO* 0 . 6?31 47606 01 
MIDDLE OF REGION 
XO. 0.4ft9?l30ME 01 
RIrHT slop OF »FR*' N 
XO. 0.54B99H31E 01 
QR * D. ij, 

QRDA* 0. 

LEFT SIDE OF RFGtO ‘ 
XO* 0,509490956 01 
MIDDLE OF REGION 
XO* 0.46620671E 01 
RIGHT Sine OF HEGl'iN 
XO* 0,53.62005* 01 
QR *0* 0. 

QRDA* J. 


J* 13 

m* ft . 19209437* ftl 

TD* 0 .12*5591 7* 01 

TD* ft,95S817niF 00 

QRADS" ft. 

QRDAS* 0. 

J. 14 

TO. 0 ■ l9t 09638F 01 
TQ* ft . 124291* 4F 01 
Tfl. ft , 909096 n 7p u 0 

QRADS* 0. 

QRDAS* 0. 
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LEFT S I DP n F RpG!0 


J" ? 0 


J» 25 


XO" 0.496*60976 01 
MIDDLE OF REGION 
XO" 0,4*i03020F 01 

RIGHT Slop or MEG I "N 

XO" 0.52490699F 01 
OR * D" U. 

QRDA" 0 . 

left Side of hfgto 
XO" 0.4R430112E 01 
H I DDL £ OF REGION 
XO" 0.43*6119*6 01 

right slop or MfcrJ 1 n 

XO" 0.5t*M54nt 01 
or*o« n, 

ORI1A" 0. 

LEFT SIDF OF RFGIO 
XO" 0.4796635*6 01 
MIDDLE OF REGION 
XO" 0.42>H971f»p 01 
RIGHT Slnp OF °FG l 'N 
XO" 0 . 50326767E 01 
GRAD" ii. 

QRDA" 0. 

LEFT SIDF OF RCGTO 
XO" 0 ,4(Sl 60071B 01 
HlnQLE OF REGION 
XO" 0.409832B6E 01 

RIGHT Slnp OF »FGl‘N 
XO" 0.49S25959E 01 
ORAD" 0 . 

OROA" 0. 

LEFT SIDE OF RFGIO 1 
XO* 0 . 451 07235F 01 
NIOOLE OF REGION 
XO" 0.39737075F 01 
rIqhT sinF of mfg^'n 

XO" 0.4R*74504fc 01 
ORAD* 0. 


TO" 0 , 19014793F 01 

TO" 0.12024574F 01 

yn« 0.B4862743F 00. 
QRADS" 0. 

QRDAS" 0. 

J« 16 

TO" 0.18927591F 01 

TO" 0 . 1164 01 "OF 01 

TO" 0.797914A3F 00 
ORADS" 0. 

QROAS" 0. 

J" 17 

TO* 0 » 1884 0 61 IF 01 

TO" 0 • 1 1274539F 01 

TO" 0.75029042F 00 
ORAOS* 0. 

QRDAS* 0. 

J" IB 

TO* D , 187614] OF 01 
TO" O.10926243F 01 

TO" 0 . 70492824C 00 
ORAOS* 0. 

ORDAS" 0. 

J" 19 

TO" 0.186B7663F 01 

TO" 0 . 1089400 OF 01 

TO" 0 , 66180676F 00 
ORAOS" 0. 


LEFT SIDF OF RPGIO 
XO" 0.44104131E 01 
MIDDLE OF REGIhn 
XO* 0.38546720F 01 
rIg h t S I df of regJ in 
XO" 0.47468689F 01 
OR AO* o. 

QRDA* ii. 

left SIDE of rfgio 

XO" 0.42?3390AE 01 
MIDDLE OF region 
XO" 0.36317123E 01 
RIGHT S I Of OF °FG I n N 
XO" 0.497B4804E 01 
ORAD" 0. 

QRDA" 0. 


LEFT S t DF OF RFGtO 
XO" 0 , 4052581 AE 01 
MIDDLE OF REGION 
XO" 0.3497060 3 F 01 
RIGHT SIDF OF MEG I 'N 
XO" 0.4495U8RE 01 
ORAD" (1, 

QRDA" 0. 

left side of rfgio 
XO" 0 . 38q 5921 OF 01 
MIDDLE OF REGION 
XO" 0.32-I86282F 01 
RIGHT SIDF OF REGION 
XO" 0.42849297F 01 
ORAD" n. 

QRda« ii. 

LEFT SIDE OF RFGIO' 
XO* 0.37*176406 01 
MIDDLE OF REGION 
XO" 0 . 30645020F 01 
RIGHT SIDF OF REGION 
XO" 0 . 41*6298 OE 01 
ORAD* (i. 

QRDA" 0. 

AND AVR 


TO* 0.18618354F 01 

YD* 0.10276649F Ol 

YD" 0.62081493F 00 
QRADS" 0. 

QRDAS" 0. 

J" 21 

YD" 0.18497644F 01 

TO" 0.96822329F 00 

TD. 0 . 54421668F 00 
ORADS" 0. 

QRDAS" 0. 

J" 22 

TO* D . 183940 7*F 01 
TO* 0 • 91366258F DO 

TD" 0.474339S9F-00 
QRADS" 0. 

QRDAS" 0. 

J" 23 

TO" 0 » 1830924 OF 01 

YD" 0.86342612F 00 

TO" 0.410360S9F-00 
QRADS" 0. 

QRDAS* 0. 

J" 24 

TD" 0.18228347F 01 

YD" O.8170036JF 00 

TD" 0. 351634176-00 
QRADS" 0. 

QRDA$* 0. 


OROA" 0. QRDAS" 0. 

4!TH BLOCkAGF rOPR*-CT 1 ON, T^E RADIATIVE F| |jV SY THE AVT 
8Y THE ABSORPTION <ETHOD. IT IS 0. 


MfTMQO IS 0. 


LEFT SIDE OF RCGIO' 

XO" 0.36A39194E 01 TO" 0.18196413E 01 
middle OF REGION 

XO" 0.29*??514E 01 TO" 0.79507545F 00 
"W * 1 ne °F t, EG I ,N 

XO" 0.4096030TF 01 TO" 0 . 32387094C-0 0 
ORAD" 0 . QRADS" 0. 

QRDA" 0 . QRDAS" 0. 


•TU/8EC-80 FT 


FIGURE 1 1 . RADIATION PROGRAM, SAMPLE PRINT OUT (Concluded) 
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Mach Number 
(FM(J)) 


Theta 
(Radians ) 


Theta 

(Degrees) 


l.'l 

1.4 
1.6 
1.8 
2.0 

2.4 

3.0 

4.0 

5.0 
10.0 


.02827 

.18325 

.28796 

.38045 

.45724 

.57417 

.68935 

.79232 

.84293 

.91972 


1.62 

10.5 

16.5 
21.8 
26.2 
32.9 

39.5 
45.4 
48.3 
52.7 


NOTE: Values at 1.4, 1.8, 5.0, and 10.0 are interpolated from 

Figure 12. Notice that because of the purpose of the table 
our values for the maximum Theta at any given Mach number are 
on the low side of the actual maximum. 


FIGURE 13. MAXIMUM VALUES OF THETA FOR SELECTED MACH NUMBERS 
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This is a typical intersection region (l.R.).The region minor axi 5 is on a line 
joining the rocket engine centers and the major axis is perpendi cular to it as shown. 



“O — 

O n3 
a> £ 
u ^ 

X 

a> a> 
Lf 

O 

> in 
J) >. 00 
C ro — 
2 o 

^ — CD 
</» fD 

E ui c 
— <u 

X +J ^ 
Q U_ 
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FIGURE 14. GEOMETRY OF INTERSECTION REGION LOCATION 





FIGURE 16. BASE GEOMETRY FOR SAMPLE PROBLEM 





CHAIM 1 TNP13X. 


1 Y U^77 1-Z. ; 1.<HiAmOO.i-,7Aii.,i. i - > A * 

♦ Block Z 

.f|1Q S* 




i53300,«^0~‘77ji3*7V A r 1 * - ^.34 ♦ ♦05794 « i.QSntt* 

Block 3 

1 • 14699 ,1.20685*1 .27162*1 *334o7*i *43165 *1.52 963, 1.68460,1*79 133 *1.90591* 

•> _ f . 1 i\ 7 k u. _ > _ > n s i a-, « y . -4 i ./ « \ * s - lli i7a4.. T9 3*/«6743.2 ♦ &22-^ 5* 

* 

2.96077,3.08055*3.19207*3.33044* 

* 

5.55 *5.676*6. 121 *6.027 *6.01* 5 . 69*3,. 674*5. 640*5*80»5. 78*5. 75*5. 723 » 5 • 7 13* 
t 'i710r Ci T^ y r t ' r ^ *69 * -5. 7w t 3 • 7 a* . s - 7 oh * 4 . 7-l£*-5..-7 18.* 6.-7 18.5. 7.1r6* 5 . 7Q 1 * 

# BT6ckTT 

5 .705 *5. 7 17 *5.72 *5 .722,3.725 *5.719 *5.685* 

* 

0.59932*0.61166 ,0.69623 *0 .65766 *0.63600 *0.60129 .0.58471 *0.56404 *0.53164* 

0 1 si Qrtfo t p t 3Q391 ,0.48 5 49 ,u. 46 9 3 0-* 0 + 43436 *0.4 33804 0.41471 ,Q . 3fl945 * Q. 37429* — 

Block 5 

* 

0.36038 ,0*34923 *0.33942*0 .3310 7. *0*32416*0 .32047,* 0*32655 *0*32794 *0 • 31689* 

n. , AiiKin.n.'^w'l7.n.>y«>7.iil')yMhi .(i. >>lnnh* 

* - - - 

‘UQ^^f0t g 71f I t0^ t lUjf A ^242 + 1 .260* 1 .292* 

Block 6 

* 

1.323,1.335*1.388 ,1.437,1.486*1.361,1.609 ,1.660 *1.710 *1.743 *1.763 *1.825* 
l .«si • ) i . i .QS‘5./-Ljnh./.nbn»>.nHh.^. 117.2.165* 

* 

1 

* Black 7 

1.90, 2^ 00 *2. 10 *2 .20,2 .30 *2.40 *2.60 *2. 80 *3. 00, 3 *20 *3.30* 

* Block 8 

1.1, 1.4*1. 6*1. 8*2. 0,2. 4*3. 0*4. 0*5. u*l 0 • 0* 

• Block 9 

• 0282 7*.l6 32 5, .287 96* .3 6043 »• 45 7^:4 * .37417*. 06 933*.792 32 , .64293 » .91972* 

CHAIN fl INPUT 

. 

♦Block 10 

6.666..63.3.45044* 

♦Block 11 

-.81404,1.57079,2.09559* 

fl^ATN 4 TNPTTT . — 


•Block 12 

X6*-6. 19*0.0* 

* Block 13 

0.0,-6.19.6.19,6.19,-6.19,0.0* 

♦Block 1A 

0.0.6.19.6.19,-6.19.-6.19.0.0* 

♦Block IS 

— 3.4*0. 0*0.0»0. 0*0 .0*0.0* 

•Block 16 

10.25*3.33, 3.3 3 . ,5.33*3.33*3.33* 

% 


A Block 17 1.0. 1.0* 


FIGURE 18 


CHAIN INPUT DATA FOR INTERSECTION REGIONS 



CHAIN 2 INPUT 
— SECOND - 


*B l ock 10 6^#06 6 » • — — 


CHAIN A INPUT 

« 

-•*8 lock— 12 
““Block 13 




1 Q . ^ . 1 Q . . 

-Ufa 1 Q . 

- i \ - r\* 

^Heek 14 
«Block 15 

. IQ aA.lQ. 1 Q__fa .1 CJ.il A*'' 



1 0 f l K f 4 - A A t '* T A * t A • 3 '* \ 

f A , ^ A , 

f 4* ^3* 

17 

■bap® X t 

1«0»1»0* 




FIGURE 18. CHAIN INPUT DATA FOR INTERSECTION REGIONS 
(Concluded) 
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APPENDIX A 


With reference to Figure 19, the form factor F from the region to 
area dA^^ is given by Reference 2 to be 


dA 


§ 


(yp - yQ dx p - (Xp - xQ dy P 

2rtr 2 * 


(33) 


Since A 2 is assumed parallel to the xz plane then dy 2 = 0 and 


V-A,- § ISP • 


(34) 


By choosing dA L at the origin 


y-L = Xi = Zi = 0. 


(35) 


Then, 


„ _ K sQ dx P 

F d Al -A £ ' 2; y 7^ 


(36) 


where 


K = y 2 . 


(37) 


Now the magnitude of the radius vector, r, is defined by 


r 2 = x 2 2 + y 2 2 + z 2 2 = x 2 2 + K 2 + z 2 2 


(38) 


83 



but 


z 2 = m x 2 + b. 


(39) 


Therefore, 


r 2 = (1 + m 2 ) x 2 2 + (2mb)x 2 + (K 2 + b 2 ) 


(40) 


and 


dA A 2 




dx- 


(1 + m 2 ) x 2 2 + (2mb)x 2 + (K 2 + b 2 ) 


(41) 


Integrating and proceeding around the trapezoid in the direction indicated 
yields 


I2 “ ^ 


2jt 


/k 2 + (b x + b 2 ) 2 ' 


tan - * 1 x 



K 2 + (b + b 2 ) ‘ 


(42) 


+g 


I 

14 2 * 


'K 2 + (b x + b 2 + b 3 ) 


2 1 



and 


"dA.-A, 


(4*) 
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